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X> ■ Abstract 

Applications of relativistic light front dynamics to computing the nucleonic 
and mesonic components of nuclear wave functions are reviewed. In this ap- 
proach the fields are quantized at equal values of r = x° + x 3 . Our method is 
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to use a Lagrangian, and its associated energy-momentum tensor T ^ to de- 



l/"") , fine the total momentum operators P^ with P + as the plus-momentum and 



P~ the r-development operator. The motivation for unusual treatment of 
nuclear physics is the desire to use wave functions, expressed in terms of plus- 
momentum variables, which are used to analyze high energy experiments such 
as deep inelastic scattering, Drell-Yan production, (e,e') and (p,p') reactions. 
This motivation is discussed, and a simple overview of light front dynamics is 
presented. Some examples of ordinary quantum mechanics are solved to show 
^ \ that the formalism is tractable. The necessary quantization is reviewed and 

applied to a series of problems: infinite nuclear matter within the mean field 
approximation; a simple static source theory; finite nuclei using the mean field 
^ . approximation; low-energy pion-nucleon scattering using a chiral Lagrangian; 

nucleon-nucleon scattering, within the one boson exchange approximation; 
and, infinite nuclear matter including the effects of two-nucleon correlations. 
Standard good results for nuclear saturation properties are obtained, with a 
possible improvement in the lowered value, 180 MeV, of the computed nu- 
clear compressibility. A complicating feature of our light front dynamics is 
that manifest rotational invariance is not used as an aid in doing calculations. 
But for each of the examples reviewed here, manifest rotational invariance 
emerges in the results of the calculations. Thus nuclear physics can be done 
in a manner in which modern nuclear dynamics is respected, boost invari- 
ance in the z-direction is preserved, and in which the rotational invariance 
so necessary to understanding the basic features of nuclei is maintained. A 
salient feature is that u>, a and tt mesons are important constituents of nuclei. 
It seems possible to find Lagrangians that yield reasonable descriptions of 
nuclear deep inelastic scattering and Drell-Yan reactions. Furthermore, the 
presence of the a and co mesons could provide a nuclear enhancement of the 
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ratio of the cross sections for longitudinally and transversely polarized virtual 
photons in accord with recent measurements by the HERMES collaboration. 
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I. INTRODUCTION AND MOTIVATION 



The need for a relativistic methodology that is broadly applicable to nuclear physics has 
never been more apparent. One of our most important sets of problems involves understand- 
ing the transition between the hadronic (baryon, meson) and the underlying (quark, gluon) 
degrees of freedom. Using a relativistic formulation of the hadronic degrees of freedom is 
necessary to avoid a misinterpretation of a kinematic effect as a signal for the transition. 
In particular, the goal of understanding future high energy and momentum transfer studies 
of nuclear targets using exclusive, nearly exclusive or inclusive processes can only be met 
through using relativistic techniques. The light front approach of Dirac |1| in which the 
time variable is taken as t + z and the spatial variables are t — z,x,y 0- ||. is one of the 
promising approaches because the momentum canonically conjugate to t — z, p + = p° + p 3 , 
is directly related to the observables. 

To be specific, consider the lepton-nucleus deep inelastic scattering observed by the 
EMC experiment This showed that there is a significant difference between the parton 

distributions of free nucleons and nucleons in a nucleus. This difference can interpreted as a 
shift in the momentum distribution of valence quarks towards smaller values of the Bjorken 
variable x. In the parton model, x is the ratio of the plus-momentum k + = k° + k 3 of a 
quark to that of the target. If one uses k + as a momentum variable, the corresponding 
canonical spatial variable is x~ = x° — x 3 and the time variable is x + = x° + x 3 = t [[| . The 
r-development operator, which plays the role of the light front Hamiltonian is P~ = P° — P 3 , 
in which P^ is the total momentum operator. In equal time dynamics, the eigenstates of the 
Hamiltonian P° are complete. In light front dynamics the eigenstates of P~ are complete: 

P ~\ n ) = Pn\ n ) (!•!) 

J = 5»<n|. (1.2) 

n 

To do calculations in this framework is to use light front formalism or light front dynamics. 
We have been attempting to derive the properties of nuclei using this light front formalism. 

Light front dynamics applies to nucleons within the nucleus as well as to partons of the 
nucleons, and this is a useful approach whenever the momentum of initial or final state 
nucleons is large compared to their mass ||. For example, this technique can be used for 
(e, e'p) and (p, 2p) reactions at sufficiently high energies. It is important to realize that the 
mere use of light-front variables for nucleons in a nucleus is not sufficient to guarantee a 
reasonable result. Shouting relativity at a problem will not make it go away. It is also 
necessary to include all the relevant features of realistic conventional nuclear dynamics. 
Combining the two aspects of relativity and realistic physics provides the technical challenge 
addressed here. 



A. An Example 

A simple example involving quasi-elastic scattering of high energy electrons from nucleons 
in nuclei shows how using the light-front approach leads to important simplifications. To 
be specific, use coordinates in which the four-momentum q of the exchanged virtual photon 
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is (v, 0, 0, —\/Q' 2 + v 2 ) ■ Here Q 2 = — q 2 , and v 2 are both very large, but Q 2 jv is finite (the 
Bjorken limit). Use the light-cone variables, so that q + « Q 2 /2u = Mx (M is the mass of 
a nucleon), q~ « 2i> — Q 2 jlv, and g~ ^> g + . 

The scattering cross section for e + A — > e' + (A — 1)/ + p, where i, f represents the 
initial, final nuclear eigenstate of P~, and p the four- momentum of the final proton, takes 
the form 

~ X)/ir/ d4 P 5 (P 2 - M 2 )5^(q + p t - Pf -p)\(pJ\ J{q) \ i) | 2 , (1.3) 

with the operator J(q) as a schematic representation of the electromagnetic current. Kine- 
matic factors and terms involving the cross section are ignored here in order to bring out 
the salient features. Performing the four- dimensional integral over p leads to the expression 

do~Y.I ((pi - Pf + qf - M 2 ) | (p = Pl -p f + qJ | J( q ) \ j) \ 2 . (1.4) 
/ p f 

Equation (|1.4j) can be evaluated most easily by using a variable 

k=pi-p f , (1.5) 

with k = (k + ,k±), the momentum that the struck nucleon had in the initial state. The 
simplification emerges from the argument of the delta function, under the Bjorken limit. 
One finds: 

(Pi ~Pf + qf -M 2 ^ -Q 2 + q-(pi - p f ) + = -Q 2 + q-k + . (1.6) 
Using this approximation, and changing the variables of integration to k gives: 

^ ~ E / {-Q 2 + <l~ k+ ) \(P = k + qJ I J(q) I •> I 2 • (1-7) 

For large values of u, one may use the single-nucleon approximation: 

(p = k + qJ\J(q)i) = (f\b k \i), (1.8) 
where bh is a nucleon destruction operator. Thus we find 

da ~ E / ^f^^ + ?^ + ) I (/ I ^ I i) I 2 • (1-9) 

The important advantage is that pj does not appear in the argument of the delta function, 
or anywhere else in Eq. ( |1.9| ). Thus the sum over intermediate states can be replaced by 
unity, using Eq. (|1.2| ). In the usual equal-time representation, the argument of the delta 
function is — Q 2 + IviEi — Ef). The energy of the final state appears, and one can not do 
the sum over states. Using light front dynamics, we re-write Eq. (|1.9|) as 
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r d 2 k dk+ 

da~J a J^5{-Q 2 + <TA;+)(zfe|,). (1.10) 

Only a ground state matrix element is required; this is the simplification we seek. Then 
integration over k + and noting that the operator b\.bk is a number operator 

da ~ d 2 k ± p(Mx,k ± ), (1.11) 

where p(Mx, k±) = (i\b\+ =Mx ^ bk+=Mx,k±\i) ls ^ e probability for a nucleon in the ground 
state to have a momentum (Mi, k±). Integration in Eq. ( |1.11| ) leads to 

a ~ J d 2 k ± p(Mx,k ± ) = f(Mx), (1.12) 

with f(Mx) as the probability for a nucleon in the ground state to have a plus momentum 
of Ms, or the nucleon distribution function. The use of light-front dynamics to compute 
nuclear wave functions should allow first-principles calculations of f(Mx). Using light- 
front dynamics incorporates the experimentally relevant kinematics from the beginning, and 
therefore is the most efficient way to compute the cross sections for nuclear deep inelastic 
scattering and nuclear quasi-elastic scattering. 



B. An Apparent Puzzle 

Since much of this work is motivated by the desire to understand nuclear deep inelastic 
scattering and related experiments, it is worthwhile to review some of the features of the 
EMC effect One key experimental result is the suppression of the structure function for 

x ~ 0.5. This means that the valence quarks of bound nucleons carry less plus-momentum 
than those of free nucleons. This may be understood by postulating that mesons carry a 
larger fraction of the plus-momentum in the nucleus than in free space PJlO|. While such a 
model explains the shift in the valence distribution, one obtains at the same time a meson 
(i.e. anti-quark) distribution in the nucleus, which is strongly enhanced compared to free 



nucleons and which should be observable in Drell-Yan experiments |TTJ. However, no such 



enhancement has been observed experimentally ||12|| , and this has been termed as a severe 
crisis for nuclear theory in Ref. 



C. More Motivation 



The use of light-front dynamics should allow us to compute the necessary nuclear meson 
distribution functions using variables which are experimentally relevant. The need for a 
computation of such functions in a manner consistent with generally known properties of 
nuclei led me to try to compute nuclear properties using light front dynamics. There are 
other motivations for using the light front formalism that have been emphasized in many 
reviews [|J. One key feature is the belief (or hope) that the vacuum of the theory is trivial 
because it can not create pairs. Another is that the theory is a Hamiltonian theory and the 
many-body techniques of equal time theory can be used here too. I also quote the review by 
Geesaman et al. J7| "In light front dynamics LFD, the particles are on mass-shell, and there 
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are no off-shell ambiguities. However, ... we have little or no experience in calculating the 
wave function of a realistic nucleus in LFD" . The aim here is to provide such wave functions. 

I also like to say: Ask not what the light front can do for nuclear physics; instead ask 
what nuclear physics can do for the light front. This is to provide a set of non-trivial four 
dimensional examples with real physics content. 

There are have been several serious efforts devoted to solving QCD using light front 
dynamics |I^rB|. This is proving be a formidable task. Much of the technical difficulties 
are associated with the vanishing gluon mass and the nearly vanishing values of the current 
quark masses. In equal time dynamics, the vacuum is known to contain both quark and 
gluon condensates, and recovering the physics of the complicated vacuum using light front 
dynamics requires a careful treatment of the constraint equations and the related zero-modes. 
This is discussed nicely in the conference proceeding [JPT| . See also the most recent review 
by Brodsky et al. [@f. 

As an example of the possible difficulties one faces in solving a gauge theory using 
light front quantization, recall the work of Casher |18] . As explained by Paston et al |20| . 



Casher found that the use of the gauge A + = is necessary. In that gauge the Feynman 
propagator has an extra power of the plus-momentum in the denominator, which strengthens 
the infrared singularities and introduces the need for a special regularization. One may " cut 
out" the infrared regime momentum space, but this procedure breaks the Lorentz invariance 
until the regularization is removed. In principle, one can preserve the gauge invariance by 
using a discrete light cone quantization [|J, but one also needs to keep the zero modes. 
This procedure typically introduces enormous complications. But not all theories are so 
complicated. In the case of Yukawa theory, which is very similar to our Lagrangians, the 



difficulty with the infrared regime as already been solved in an convincing manner Ref. |T9 

The essential simplification encountered here is in our use of hadronic Lagrangians, which 
are not gauge theories. For example, consider that, in the theories we consider, all of the 
particles are massive. No one believes that the vacuum is a condensate of nucleon-anti- 
nucleon pairs. The vacuum really should be trivial. The implication is that the complications 
associated with zero modes can legitimately be ignored. Furthermore, the meson masses are 
much larger (> 135 MeV ) than typical energy scales (~ 10 MeV) of nuclear physics. This 
means that truncation of the Fock space to terms involving only 0,1 or 2 mesons should be a 
reliable approximation. Thus many of the technical difficulties of using light front dynamics 
are not present for hadronic Lagrangians. There is no reason to prevent light front dynamics 
from becoming a standard tool of relativistic nuclear physics. 



D. Outline and Scope 

The general goal is to provide a series of examples showing that light front dynamics can 
be used for high energy realistic and relativistic nuclear physics. The goal is to present the 
physics in as simplified form as possible. 

We shall begin with a very brief overview of light front dynamics. Then some simple 
examples [21], requiring little technology are solved. The formal procedures of light front 
quantization of a hadronic Lagrangian C are discussed next. Our first application of this 



quantization is a study of infinite nuclear matter within the mean field approximation [22 
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The distribution functions f(y) for nucleons and mesons are computed. The results are 
startling -the vector mesons are found to carry a substantial fraction of the nuclear plus- 
momentum. This finding is analyzed using a simple static model p3] in the following section. 
Then finite nuclei p4| are studied using the mean field approximation. Here one confronts 
the problem that the use of x~ = t — z as a spatial variable violates manifest rotational 
invariance. We found p4| that rotational invariance re-emerges in the results, after one 
does the appropriate dynamical calculation. The desire for a realistic approach means that 
nucleon-nucleon correlations must be included and one must do physics beyond the mean 
field approximation. Chiral symmetry and a simple example involving pion-nucleon scat- 
tering is then discussed |[22|| . This is necessary preparation for the development of a new 
light-front one boson exchange potential |25 . The influence of nucleon-nucleon correlations 



on the properties of nuclear matter is studied by making the necessary light front calcula- 
tions. Applications are to compute the nuclear pionic content and to nuclear deep inelastic 
scattering and Drell-Yan processes. This study indicates that the apparent puzzle is not 
a real puzzle. Then some plans for future research are discussed. The impact of the calcu- 
lations discussed here is that light front dynamics can be used for high energy realistic and 
relativistic nuclear physics. Reasonable results for standard properties are obtained. There 
is a strong desire for an experimental measurement of an observable that has an importance 
more readily appreciated and that is more easily computable using light front dynamics. As 
discussed, in Sect. 11, the discovery of the HERMES effect |26| that the ratio R = o L jo T is 
significantly enhanced in nuclei may provide such an observable. 

The emphasis of this review is on the influence of light front quantization on nuclear 
physics. Hadronic Lagrangians are used. There are several excellent reviews of light front 
techniques ||J|, but these are mainly devoted to applications in particle physics fl4Hl6l . I 



am concerned here mainly with heavy nuclei. Applications of light front techniques to nuclear 
reactions have been pioneered by Frankfurt and Strikman |2|||. The present work reviews 
the first attempts, using realistic inputs, to apply light front techniques to the structure 
of heavy nuclei. This is the salient feature of the present review. However, it is worth 
mentioning the differences between the present and previous uses of light front techniques. 



Apart from Ref. |p7| , almost all of the early work P,p|,p8|-[3T[] is devoted to the two-nucleon 
system, which is of tangential importance here. But this is only one difference. There are 
also differences in philosophy. Our approach is to choose a Lagrangian, and then determine 
its consequences for experimental observables. If the consequences are not palatable, we 
change the Lagrangian. The approach of |2%|-|5U|] relies on symmetries which are used to 
write down the allowed form of operators. This is not sufficient to determine uniquely, for 
example, all of the components of the electromagnetic current operator J M . The approaches 
of [[3(],[31]] try to restore manifest rotational invariance by introducing a general direction n to 
replace the three-direction. We expect, and show, that the final calculations of observables 
do respect this invariance. Our approach is closest to that of Refs. which is based on 
Feynman diagrams. The specification of a Lagrangian, should allow us to go further because 
all of the relevant diagrams can be defined. 



8 



II. WHAT IS LIGHT FRONT DYNAMICS? 



This is a relativistic treatment of dynamics in which the fields are quantized at a fixed 
"time'V = t + z = x° + x 3 = x + . This means that the independent spatial variable must be 
x~ = t — z so that the canonical momentum is p° + p 3 = p + . The remainder of the spatial 
variables are given by x±,p±. 

The consequence of using r as a "time" variable is that the canonical energy is p~ = 
p° — p 3 . Our notation for four- vectors is to use: 

A ± = A°±A 3 , (2.1) 

with 

A ■ B = A^B^ = - [A + B- + A~B + ) -A ± -B ± . (2.2) 

The key reason for using such unusual coordinates is phenomenological. For a parti- 
cle with v ~ ce3, the quantity p + is BIG. Thus experiments tend to measure quantities 
associated with p + . 

Another important feature is the relativistic dispersion relation p^p^ = m 2 , which in 
light front dynamics takes the form: 

P- = Pl ^~. (2.3) 
p+ 

Thus one has a form of relativistic kinematics that avoids using a square root. Eq. ( |2.3|) is 
perhaps the most important equation in light front physics. Its consequences appear again 
and again in this review. If a reader, new to the light front, decides to learn only one of the 
many equations of this review-this should be the equation. 

The main formal consequence of using light front dynamics is that the minus component 
of the total momentum, P~ , is used as a Hamiltonian operator, and the plus component 
P + is used as a momentum operator. These are obtained by using the energy-momentum 
tensor 

dC 

T^ = -g^£ + ——d» ( p i , (2.4) 

in which the degrees of freedom are labeled as 0j. The term T +At is the density for the 
operator P M , with 

P» = - J d 2 x ± dx-T + ». (2.5) 

The use of a Lagrangian to provide these operators is the salient feature of our approach. 
The procedures to obtain these operators for the different Lagrangians are discussed below. 
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III. SIMPLE PROBLEMS 



We introduce light front techniques by solving the wave equations for the three dimen- 
sional (scalar) harmonic oscillator, and the three dimensional Hulthen potential. We con- 
sider a situation in which a single nucleon moves under the influence of an external potential. 
This external potential arises from interactions with residual (A — l)-particle nucleus. The 
physics of the entire A-nucleon system is Lorentz invariant. In the present example the 
external potential is most simply expressed in the rest frame of the nucleus, for which it 
independent of time. For simplicity, we neglect the effects of spin and use the Klein-Gordon 
equation. 



A. Three Dimensional Harmonic Oscillator 

In ordinary coordinates, the Klein-Gordon equation for a 3 dimensional relativistic os- 
cillator 



p + m + kx (f) (3-1) 



closely resembles the Schroedinger equation for a non-relativistic harmonic oscillator, and 
the energy eigenvalues are 

E% = m 2 + u(n+^, (3.2) 

where us 2 = 4k. 

The light front version of the Helmholtz equation is obtained by changing variables to 
x = x° ± x 3 , with 

Then 

(<T<9 + - Vi + m 2 + k{x\ + (x + - x-f/A)) = 0. (3.4) 

A potential that is static in the equal time formulation corresponds to a LF-"time", i.e. x + , 
dependent potential in light-front coordinates. A static source in a rest-frame corresponds 
to a uniformly moving source on the light-front because the line z = Zq corresponds to 
x ± = t ± Zq and Ax~/Ax + = 1. Thus the time dependence is only due to a uniform 
translation, and it should be easy to transform Eq. ( |3.4j) into a form which contains a 
potential that does not depend on x + . For this purpose, we consider the equation of motion 
satisfied by fields obtained by an x + (LF-time) dependent translation 

<f)(xj_,x~ ,x + ) = e~ ix+p+/2 x{x±, x~,x + ) (3.5) 
Using P + = —i2-7^pr, we find 
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D ix+P+/2 



f 



I T — X 



-P+/2 



f 



-X 



e ix+P+/2Q~ e -ix+P+/2 _ Q- _ Q+ 

The use of Eqs. ( [Ol) and (|3.6| ) in Eq. (|3.4j ) then leads to the result 



(<T - d + )d + - Vi + m 2 + k{x\ + (^-) 2 ) J X = 0. 



(3.6) 



(3.7) 



Note that the potential does not depend on x + and is instead a function of x\ + (^r) 2 . 
The term (^-) 2 can be thought of as playing the role of z 2 , so that the potential is actually 
spherically symmetric. Indeed, if we simply set x + = 0, then x° = —x 3 and x~ = — 2x 3 = 
-2z. 

It is instructive to write the result (|3.7|) in an operator form: 



id x 



-Vi + m 2 1 
id + H 7j- h 



A- 



(3.8) 



where V = k(x^ + (^-) ). This looks like an operator version of Eq. ( |2.3|) , but including 
an interaction and an additional p + operator which accounts for the recoil of the A — 1- 
nucleon residual nucleus, see Sect. VII. The factors ^= are boson normalization factors; 
see Sect. IV. 

The potential appearing in Eq. 
of the form: 



X(x) = e 

Using this and completing the square, leads to the result 



does not depend on x + so that there will be solutions 

ip - x+/2 Xn(x~,x ± ). (3.9) 



X 



-{d + + ivljlf - Vl + m 2 + k{x\ + (— ) 2 ) X n = (P"/2)V 



(3.10) 



One converts the operator (d + + ip n /2) 2 to (<9 + ) 2 using yet another transformation: 



X{x ,x± 



~) _ e -ip n x /4 



tC ^ tC _|_ J 



to find 



-(d H 



V i + m 2 + k 



F = (p;/2) 2 F. 



(3.11) 



(3.12) 



This is the same form as the equation in the equal time coordinates, and p n /2 takes on the 
values of E„: 



(p-/2) 2 = m 2 + cu(n + 3/2). 



(3.13) 



The light front version gives the same results for the eigenvalues as the ET development. 
So one might wonder why we are doing the light front at all. The point is that we are able to 
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compute the light front wave functions that depend on x , or in momentum space depend 



on p + . The wave function of the ground state of Eq. (3.4) is given by 



e ^-/4jV exp l--^(x 2 ± + ?—)]. (3.14) 



The number density no(p±,p + ) is defined as the square of the momentum space version 
of xo- This quantity is accessible in high energy proton and electron nuclear quasi-elastic 
reactions. It is useful to define the light front variable 

«^ + /(y)- (3-15) 

For a shell-model nuclear target p$ is a definite fraction of the total minus-momentum. In 
the rest frame, the minus-momentum is the same as the plus- momentum. Thus the variable 
a is a ratio of plus-momenta and uq(p±, a) is independent of frame. 

The Fourier transform of Eq. Q3.14|) is obtained by multiplying by e l P±- x ± e ~ l 2P +x and 
doing the integral over all x± and x _ /2. Then one determines that 

_ _pi -( P ~( Q -i)) 2 

no(p±,a)=N e 4 ^ . (3.16) 

Note that one finds the same p± distribution for each value of the variable a. This is not a 
general feature of light-front wave functions, as we show in the next sub-section. 

In an exact calculation, the wave functions and number density should vanish for values 
of a plus- momentum fraction variable that are not between and 1. This is referred to as 
having proper support or satisfying the spectrum condition. This property is not evident 
here because the system consists of one particle that interacts with an external source which 
supplies the external potential V and is also responsible for the term id + on the right hand 
side of Eq. ( |3.8| ). To better understand the support and the origin of the term id + , consider 
a two-particle problem in which a light particle (nucleon) of mass m interacts with a heavy 
particle (residual nucleus) of mass M. The total four-momentum of the system is denoted 
as PP. For a nucleus with A nucleons, M « (A — l)m and M 3> m. This is a weak binding 
limit (binding energy per particle much less than the mass) appropriate for nuclear physics. 
The light particle is labeled as 1 and the heavy particle as 2. The plus- momentum fraction 
carried by the light particle is denoted as (3 with 

= fk (3-17) 



so that 

P Po 



a 2P^ 



(3.18) 



Using the given masses m, M it will always be true that /3 <C 1, and that a ~ 1. The relative 
momentum is given by 

p± = (1 - P)pi± - Pp*±, (3.19) 
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with 



P± « Pi± (3.20) 

The use of these variables allows leads to a simple expression for the kinetic energy 

_ _ p\ m 2 M 2 

Pl +P2 ~ (3(1 - P)P+ + + (1-/3)P + ( J 

p 2 , +m 2 M 2 . 

^^^+P + + ~^, 3.22 

in which the weak binding approximation j^ 2 — > 1 is used. The right hand side of Eq. ( |3.22| ) 
is the kinetic part of the right hand side of Eq. (|3.8|) . Thus the term id + = p + represents 
the momentum of the heavy particle which supplies the external potential. The use of the 
kinetic energy operator of Eq. ( ft.21|) instead of that of Eq. (|3.8|) would lead to obtaining 
exact support properties. 

The approximation does not cause much trouble. For large values of the particles mass 
m, (m ^> K 1 / 4 ) (which corresponds to the situation of nuclear physics in which the product 
of the nucleon mass and the nuclear radius is a very large number) the value of a must 
always be close to unity. In nuclear physics ^ is of the order of (5-Ra/Ful) 2 which taking 
Ra = 4 Fm, yields e~ 400 at x — 0. In that case, there is no problem with the support. 

To provide an example, we consider the electromagnetic form factor of the ground state 
of a heavy nucleus. This has been measured, for many nuclei, in elastic electron scattering, 
and the sizes of nuclei have been determined as one of the classic achievements of nuclear 
physics. Interest in this topic has been revised because of a recent proposal to Jefferson 
Laboratory to use parity-violating electron scattering to measure the neutron radius [32 . 
Determining this to high precision is needed, and can be obtained provided one knows the 
proton distribution. Therefore it is useful to examine the influence of presumably small 
effects such as relativistic corrections. One works using a reference frame, the Drell-Yan 
frame [|33] in which the plus component of the four vector of the virtual photon vanishes, 
so that Q 2 = —q 2 = q\. Now consider the various momenta involved. Suppose initially, 
P± = 0. If the light particle absorbs the virtual photon the final momenta are p± + q± for 
particle 1 and — p± for particle 2. The relative momentum is then p± + (1 — (3)q± In this 
case the form factor F(Q 2 ) (matrix element of the plus component of the electromagnetic 
current operator) is given by 

F(Q 2 ) = J d 2 p ± d(3P X o(P,P±)Xo(P,P± + (l-P)q±), (3.23) 

in which the influence of relativity appears in the integral over f3 and the factor f3. For the 
harmonic oscillator ground state we find: 



P± , Q 2 (l-/3) 2 \ -2 



F(Q 2 )=N j d 2 P± daae~^ ^ ) e -^ (a - 1)_ , (3.21) 

where 
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Po 



m 2 + Zy/~K. 



(3.25) 



Our purpose here is the study of nuclear physics, so we are interested in the non- 
relativistic limit and the corrections to it. To this end, we define a variable p z using 



a 



1 + P -. 

m 



(3.26) 



The non-relativistic limit of (|3.24 ) is obtained by letting m approach infinity, but with 
m/M » l/(A- 1), and 



m 
M 



(3.27) 



Then ^- = m 2 and we find 



NR 



(Q 2 ) = N NR ^l e -^ + ^^ 

J m 



The mean square radius —6 



dF(Q 2 ) 
dQ 2 



P 2 



— g 4\/re 

Q2=o, is given by 
3(1 -m/M) 2 



(3.28) 



2 A 

The experimental value of the nuclear radius is given approximately by 

i? exp = 1.1 A 1/3 fm. 



(3.29) 



(3.30) 



The leading correction to this comes from the product of the factor — in the integrand 
of Eq. ( [PI with a ff which comes from a term in the exponent. This is correction is 
of order p 2 z /mM ~ ^/J^/mM. We define a semi-relativistic limit SR via the use of ( |3.24| ) 
and keeping the leading correction terms. Performing the straightforward evaluations, using 
Rnr = R exp , leads to the result: 



P2 p2 
r _ f^SR ^NR 
— d2 

NR 



mM' 



(3.31) 



or using Eq. ( |3.30D 



5 



0.055 



AVZ(A- 1)' 

This corresponds to very small (8 x 10~ 6 ) effects for large nuclei A ~ 200. 



(3.32) 
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B. Light Front Hulthen Wave Function 



Any static potential of the form V(x 2 = x\ + z 2 ) can be solved on the light front. The 
transformation ( |3.5| ) corresponds to including the p + term in the x + development operator 
and a simple prescription of replacing z by —x~/2 in V. We present here the solution for 
the Hulthen potential. This allows us to demonstrate an interesting contrast between the 
implications of different forms of potentials. In the equal time formulation we have the wave 
equation: 



E 2 



in which V H is chosen so that the lowest energy solution is 

0(r) = N(e~ 
where b > a. The eigenenergy is given by 



-ar g— br\ 



E = \f' 



m z — GT 



(3.33) 



(3.34) 



(3.35) 



The light front version of the wave equation is obtained in the same manner as Eq. ( [3.12D , 
with the result 



-X[X ,2_L, 



-2i- 







dx~ 



+ pl + V H [xl + 



x~ 



X{x ,x±). 



(3.36) 



The lowest value of p n /2 is clearly the same as E of Eq. ( |3.35| ), and the ground state wave 
function is given by 



,fi) = e ip " x ~ /4 N^ 



exp 



-a\ Xi + 



exp 



-b\ x 2 < + 



x~ 



(3.37) 



The momentum distribution (p±,p + ) obtained here provides an interesting contrast with 
that of the harmonic oscillator ( |3.16| ). We find 



(a 2 +p 2 ± + (p-/2) 2 (a - I) 2 ) 2 (b 2 +p 2 ± + (p-/2) 2 (a - 1) 2 Y 



(3.38) 



There is a different p± distribution for each value of a, with broader functions of p± obtained 
for smaller values of a; see Fig. |l]. The results in the figure are obtained using m=.94 GeV, 
£'=.932 GeV and b = 5a. An experimental hint of such a behavior has been found recently 
31. 



IV. LIGHT FRONT QUANTIZATION 

I discuss the basic aspects in an informal way as possible. For more details, see the 
reviews and the references. It's easiest to start by considering one free field at a time. These 
are the scalar meson 0, the Dirac fermion ip and the massive vector meson V^. 
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p,,GeV 

FIG. 1. Light front momentum distribution as a function of a and p±. 



A. Free Scalar field 



Consider the Lagrangian 



1 



2 a3\ 



The notation is such that d ± = d° ± d 3 = 2^^. The Euler-Lagrange equation leads to the 
wave equation 



id 



-Vi + m 2 s 



The most general solution is a superposition of plane waves: 



x 



d 2 k ± dk + 6(& 



a(k)e~ ik - x + a\k)e 



Ak-x 



where k-x = ^(k~x + + k + x^) — k± • x± with k~ - ,J —- 
restricts k + to positive values. Note that 



(4.2) 



(4.3) 



k+ ' 



id+e -ik-x = k + e ~ik-x_ 



and k = (k + , k±). The 6 function 



(4.4) 



The value of x + that appears in Eq. (|4.3|) can be set to zero, but only after taking necessary 
derivatives. 

Deriving the equal x + commutation relations for the fields is a somewhat obscure pro- 



cedure |31>H37L but the result can be stated in terms of familiar commutation relations: 



[a(k),a\k')} = 5(k ± - k' ± )5(k + - k' + ) 



(4.5) 
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with [a(fe),o(Jfe')] = 0. 

It is interesting || to consider the commutation relations between fields at the same 
value of light front time: 

[<f>(x), <P(y)} x+=y+ = ^e(x- - y-)6®(x ± - y±), (4.6) 

where e(x) = 8(x) — 9(—x). This commutator, evaluated for x° = y°, vanishes because 
the separation between x and y would necessarily be space-like. However, for x + = y + the 
separation is — (x± — y±) 2 which vanishes for x± = y±. 

The next step is compute the Hamiltonian P~ for this system. The conserved energy- 
momentum tensor is given in terms of the Lagrangian: 

T *^-^ + oiM (T4 - (47) 

This brings us to the question of what is g^ u l This is straightforward, although the results 
(viewed for the first time) can be surprising: 

9 +v = 9° v + 9 3u , (4.8) 

which implies 

g ++ = g 00 + g 03 + g 30 + g 33 = 1 + + — 1 = 
9 ij = = 1,2, J = 1,2); g + ~ = g' + = 2. (4.9) 

Then one finds that 

= Vj_0 • V_l0 + m 2 s (f) 2 . (4.10) 
In general, the term T + ~ is the density for the operator P~: 

P~ = -J d 2 x ± dx~T + -. (4.11) 



The use of the field expansion Q4.3|) , along with normal ordering and integration leads to the 
result: 

7 2 2 

p- = J d 2 k L dk + 6(k + )a\k)a(k) ± ^ + ms - (4-12) 

One defines a vacuum state | 0) such that a(p) | 0) = 0. Then the creation operators acting 
on the vacuum give the usual single particle states: 



^ t b)|0) = ^-^at(p)|()). (1.13) 



V 

The momentum operator P + is constructed by integrating T ++ : 



p+ = J d 2 k ± dk + 6{k + )a\k)a{k)k + . (4.14) 
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It is interesting to consider how the results for P~ are changed if interactions are included. 
If we add an interaction term of the form \(pV(j) one would get an equation for P~ of the 
form of Eq. ( |3.8| ), but without the term id + on the right hand side. 

Another consideration involves the vacuum. Suppose we take the Lagrangian 

C = ^(d^dy - m 2 s <p 2 ) + A0 4 . (4.15) 

The operator creates or destroys a particles of plus- momenta k + > 0. Thus possible terms 
in which A0 4 term converts the vacuum | 0) into a four particle state vanishes by virtue 
of the conservation of plus-momentum. The vacuum of p + = can not be connected to 
four particles, each having a positive k + . This vanishing simplifies Hamiltonian (x + -ordered 
perturbation) calculations. However, one must be careful about physics involving k + = 0. 
For example, Burkardt constructs non-vanishing perturbative corrections to the self 
energy of the scalar meson which do involve the spontaneous creation of four scalar mesons 



from the vacuum. For a detailed discussions, see Ref. [17|. 



B. Free Dirac Field 

Consider the Lagrangian 

A^^-^-M)^, (4.16) 

and its equation of motion: 

- Af) V = 0. (4.17) 

A fermion has spin 1/2, so there can only be two independent degrees of freedom. The 
standard Dirac spinor has four components, so two of these must represent dependent degrees 
of freedom. In the light front formalism one separates the independent and dependent 
degrees of freedom by using projection operators: A± = |7°7 ± . Then the independent field 
is ip + = A + ip and the dependent one is ?/>_ = A_-0 
The Dirac equation fl4.17|) is re-written as 

2 7 + <T + -7~<9 + + n ± ■ Vj_-Af)v = 0. (4.18) 



2 ' 2 



Equations for ijj± can be obtained by multiplying Eq. ( |4.18|) on the left by A 



id~i)+ = (a ± ■ — + (3M)i)_ 
i 

id + 4>_ = ( a± - — + (4.19) 

t 



so that the equation of motion of ip + becomes 



id~^+ = ( a± -^r L + (3M)-^(cx ± ■ + /3M)V+. (4.20) 
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One can make the field expansion and determine the momenta in a manner similar to 
the previous section. The key results are 



1 / Vi \ 

V • -=- + PM Uh 



= E / d 2 P± dp + 6(p + ) 



id 

pi + m 2 



6 t (p,A)fe(p,A) + ci t (p,A)ci(p,A) 



(4.21) 
(4.22) 



where b(p, A), A) are nucleon and anti-nucleon destruction operators. 



C. Free Vector Meson 



The formalism for massive vector mesons was worked out by Soper |39| and later by Yan 



40| using a different formulation. I generally follow Yan's approach. The formalism in the 



references is lengthy and detailed in the references, so I try to state the minimum. There 
are three independent degrees of freedom, even though the Lagrangian depends on and 
ytxv = Qvyv _ guyfi Thege are c h osen to be V + and V +i . The other terms V~,V\ V~ l and 

V l i can be written in terms of V + and V +l . 

The expression for the vector meson field operator is 

= I ^wf e ^ <») K fc ' + at ( fc ' > ( 4 - 23 ) 

J (27r) 3 / 2 V2fc+ w=1)3 J 
where the polarization vectors are the usual ones: 

fc%(fc, u) = 0, e"(fc, u;)e M (fc, u/) = -<W, 

£ e"(fc, W )e"(fc, w) = -(<T - =4-). (4-24) 

Once again the four momenta are on-shell, with k~ = ± k+ - . The commutation relations 
are 

[a{k,u),a\k',u')] = 5^5^ + \k - k'), (4.25) 
with [a(k,u),a(k',uj')] = 0, and lead to commutation relations amongst the field operators 



that are the same as in Ref. 40 



D. We need a Lagrangian no matter how bad 

It seems to me that one can not do complete dynamical calculations using the light 
front formalism without specifying some Lagrangian. One starts [^] with £ and derives 
field equations. These are used to express the dependent degrees of freedom in terms of 
independent ones. One also uses C to derive T^ v (as a function of independent degrees of 
freedom) which is used to obtain the total momentum operators P . It is P~ that acts as 
a Hamiltonian operator in the light front a: + -ordered perturbation theory. 
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We start with a Lagrangian containing scalar and vector mesons and nucleons if)'. This is 
the minimal Lagrangian for obtaining a caricature of nuclear physics because the exchange 
of scalar mesons provides a medium range attraction which can bind the nucleons and the 
exchange of vector mesons provides the short-range repulsion which prevents a collapse. It 
is also useful, for phenomenological purposes, to include scalar meson self-coupling terms. 
Thus we take 



+ ^ (7 M (^ % ~9v V„) -M- g a ^j i>\ (4.26) 
with the effects of other mesons included elsewhere and below. The equations of motion are 
d,V^ + rr%V = g^Yip' (4.27) 
d^<f> + m 2 s <p + ~</> 2 + = (4.28) 

(id- - g v V-)if>' + = (a x ■ ( P± - g v V ± ) + (3(M + g s (f>))if)'_ (4.29) 
{id + - g v V+)iP'_ = (a_L • ( P± - g v V ± ) + (3(M + g s <f>))if}' + . (4.30) 

The presence of the interaction term V + on the left-hand side of Eq. ( 4.30Q presents a 
problem because one can not easily solve for ip_ in terms of if> + . This difficulty is handled 
by using the Soper-Yan transformation: 

if)' = e- i9 ° A{x) ip, d + A = V + . (4.31) 

Using this in Eqs. ( f4.29| )-( f4.30| ) leads to the more usable form 

{id~ - g v V~)ifj + = (a ± ■ (p ± - g v V ± ) + p(M + g s 4>))^. 

id + if>„ = (a ± ■ (p± - g v V±) + (3(M + g s (f>))^ + . (4.32) 

The cost of the transformation is that one gets new terms resulting from taking derivatives 
of A(x). One uses with 

V» = V» - 7^d>*V + , (4.33) 

and enters in the nucleon field equations, but V^ 1 enters in the meson field equations. 
We also need the eigenmode expansion for V^. This is given by 

J (2tt) 3 / 2 V2k+ w=1)3 J 
where, using Eqs. (|4.33|) and ( |4.23|) , the polarization vectors e M (fc,c<j) are 

e"(fc, u) = e"(fc, u) - ^e + (fc, w). (4.35) 

Note that 
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0J=1,3 

Then we may construct the total four-momentum operator from 

= - J dx-d 2 x ± T + ^(x + = 0,x~,x ± ), 



(4.36) 



with 



dC 



^ <9(<9 M r ) 



d»<j> r , 



(4.37) 



(4.38) 



in which the degrees of freedom are labeled by <p r . We need T ++ and T + , which are given 
by 



d + (j)d + <j) + V ik V ik + m 2 v V + V + + 2^\id + %l) 



(4.39) 



and 



A 



T + - = V ± ■ V ± + mj0 2 + 2 (^-0 3 + -0 4 



l(y + -) 2 + ~V kl V kl + m 2 v V k V k 



(4.40) 



This form is still not useful for calculations because the constrained field ip- contains 
interactions. We follow Refs. in expressing ?/>_ as a sum of terms, one £_ whose 

relation with is free of interactions, the other r/_ containing the interactions. That is, 
rewrite the second of Eq. ( f4.32| ) as 



v- 



id+ 
1 

id+ l 



{a± ■ p± + /3M)^+ 
a± -g v V± + pg s (j))ip + . 



Furthermore, define = i/; + (x), so that 

ip(x) — £(x) + rj-(x) 



(4.41) 



(4.42) 



where = + £+(x). This separates the dependent and independent parts of ip- 

One needs to make a similar treatment for the vector meson fields. The operator V + ~ 
is determined by 



V 



g v J + - m 2 u V + - diV 



(4.43) 



Part of this operator is determined by a constraint equation, because the independent vari- 
ables are V + and V l+ . To see this examine Eq ( 4.43Q , and make a definition 



v- 
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(4.44) 
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where 



d+ J ■ 



(4.45) 



The sum of the last term of Eq' ( |4.4C ) and the terms involving u + ~ is the interaction density. 
Then one may use Eqs. ( f4.40|) , (|4.42 ), and ( |4.44| ) to rewrite the P~ as a sum of different 
terms, with 



P 



ON 



d 2 x ± dx' £ (7 ± • p± + M) f , 



(4.46) 



and the interactions 



Pj = Vi + V 2 + t> 3 , 



(4.47) 



with 



J d 2 x±dx £ (g v j ■ V + g s <f) £, 



(4.48) 



V 2 



d 2 x±dx £ (-flf«7 • V + g s (f) (-^7 • v + 9s4>) € 



(4.49) 



and 



v 3 



9v 



J d 2 x±dx J dy x e(x - y x ) Z+(y x , x ± )<y + £ + (y x , x±) 
x / dy 2 e(x~ - y 2 )^ + (y 2 ,x ± )~f + £ + (y 2 ,x ± ), 



(4.50) 



where e(x) = 8(x) — 8(—x). The term v x accounts for the emission or absorption of a single 
vector or scalar meson. The term v 2 includes contact terms in which there is propagation 
of an instantaneous fermion. The term t> 3 accounts for the propagation of an instantaneous 
vector meson. 

Our variational procedure will involve the independent fields so we need to express 
the interactions P$ N and v X}2 in terms of £+. A bit of Dirac algebra shows that 



P 



N 



P QN + Vl+V 2 

d 2 x ± ^fzl[2g,v- 



+ (ol ■ (rl - g v V±) + P(M + g s 



id+ 



oc± • (p± - g v V ± ) + (3{M + g^)) £+. (4.51) 



It is worthwhile to define the contributions to P arising from the mesonic terms as P^ 
and P^. Then one may use Eqs. ( f4.40|) and ( [4.39|) along with the field expansions to obtain 



1 f 

P~ = - J d 2 x ± dx~ (v_l</> ■ V_l0 + m 2 s (, 

k 2 , + m 2 



d 2 k ± dk + 6(k + )a\k)a(k)' 



k+ 



(4.52) 
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p+ = J d 2 k ± dk + 6(k + )a\k)a(k)k + , (4.53) 



7 2 2 

P~ = Y [ d 2 k ± dk + 6(k + )^^aUk,u)a(k,Lu) +v 3 , (4.54) 
. ~? o J k + 



u=l,3 

and 



P+= I d 2 k ± dk + 6(k + )k+a\k,uj)a{k,uj). (4.55) 

, 1 Q <J 



w=l,3 



The term V3 is the vector-meson instantaneous term, and we include it together with the 
purely mesonic contribution to P~ because it is canceled by part of that contribution. 
Thus, our result for the total minus-momentum operator is 

P- = P N +P- + P~, (4.56) 

and for the plus-momentum 

P+ = P+ + P+ + P^ (4.57) 

where from Eq. ( |4.39| ) 

Pn = \J d 2 x ± dx~ 2^ + id + ^+. (4.58) 



V. INFINITE NUCLEAR MATTER MEAN FIELD THEORY 

The aim of our approach is to do the realistic, relativistic physics of large nuclei using 
light front dynamics. Since this is a difficult project, it is very worthwhile to consider first 
the simplest possible example. This is to consider the nucleus of infinite size. One considers 
a limit in which the radius and nucleon number are each taken to be infinity, but the baryon 
density, or number of nucleons per unit volume, is taken as finite. Furthermore, a useful first 
step is to consider the mean field approximation to to the dynamics in which nucleons are 
treated as moving under the influence of self-consistent potentials. Field-theoretic versions of 
relativistic mean field theory were pioneered by Walecka and exploited heavily by his students 
and many other workers |[42||. Until our work, the only light front version of relativistic field 



theory applied to infinite nuclear matter was that of Glazek and Shakin |E| who treated 
the nucleus as bound under the influence of scalar mesons. 

The philosophy |42| behind field-theoretic mean field theory is that the nucleonic densities 



which are mesonic sources are large enough to generate a large number of mesons to enable 
a classical treatment (replacing the operators which represent the sources by expectation 
values). In infinite nuclear matter, the volume is taken as infinity so that all positions 
are equivalent. Then, the nucleon mode functions are plane waves and the nuclear matter 
ground state is assumed to be a normal Fermi gas, of Fermi momentum kp and of large 
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volume Q. Under these conditions one finds solutions of Eqs. ( |4.27f ) and ( |4.28| ) in which the 
meson fields are constants: 



0=-iL<^) = _iLp s (5.1) 
Trig mj. 

where the expectation values refer to a ground state expectation value, and Pb = 2A^/37r 2 . 
This result that is a constant. The notion that there is no special direction in space is 
used. Eq. ( |4.33| ), tells us that the only non- vanishing component of V is V~ = V°. Since the 
potentials entering the light-front Dirac equation ( [4.32 ) are constant, the nucleon modes are 



plane waves ip ~ e * fc - x ; anc l the many-body system is a Fermi gas. The solutions of Eq. ( |4.32| ) 
are 

^ + = g ^ h + k l±±iL±Mt^ (5 .3, 



Solving the equations (|5.1|),fl5~2"|) and ( |5.3| ) yields a self-consistent solution. The light front 



eigenenergy [id~ = k~) is the sum of a kinetic energy term in which the mass is shifted by 
the presence of the scalar field, and an energy arising from the vector field. Comparing this 
equation with the one for free nucleons shows that the nucleons have a mass M + g s cf) and 
move in plane wave states. The nucleon field operator is constructed using the solutions 
of Eq. ( |5.3| ) as the plane wave basis states. This means that the nuclear matter ground 
state, defined by operators that create and destroy baryons in eigenstates of Eq. ( |5.3|) , is 
the correct wave function and that Equations ( |5.1| ), ( |5.2j ) and ( |5.3| ) represent the solution 
of the approximate field equations, and the diagonalization of the Hamiltonian. 

A. Nuclear Momentum Content 

The expectation value of T +fl is used to obtain the total momentum: 

P» = Iy d 2 x ± dx-(T +l *). (5.4) 

The expectation value is constant so that the volume ft = \ j d 2 x±dx~ will enter as a factor. 
A straightforward evaluation leads to the results 

^ - m y + 2 + ^) + ^ I w ( 5 „ 

^ = mliV-f + j^y 3 J p d 2 k ± dk + k\ (5.6) 
To proceed further one needs to define the Fermi surface F. The use of a transformation 



k + = V (M + g s <pf + k 2 + k 6 = E(k) + k 6 (5.7) 
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to define a new variable k 3 enables one to simplify the integrals. One replaces the integral 
over k + by one over k 3 , including the Jacobian factor 

~W - ~E- (5 ' 8) 

Then one computes the nuclear energy E as the average of P + and P~: E = (P~ + P + ). 
The results are the very same expressions as in the original Walecka models. This provides 
a useful check on the algebra, which is important because this model has been solved in a 
manifestly covariant fashion. 

There is a potential problem: for nuclear matter in its rest frame we need to have 
P + = P~ = Ma- If one looks at the expressions for P this result does not seem likely. 
However, the value of the Fermi momentum has not yet been determined. There is one more 
condition to be satisfied: 



( d(E/A) \ 
V dk F ( 



0. (5.9) 



Satisfying this equation determines k F and for the value so obtained the values of P + and 
P turn are the same. 

Thus we see that our light front procedure reproduces standard results for energy and 
density. We discuss two sets of results. The first involves the original Walecka model, in 
which A = n = 0. Then we may use the parameters of Chin and Walecka |44| g 2 M 2 /m 2 = 
195.9 and g 2 M 2 /m 2 = 267.1 to obtain first numerical results. In this case, k F = 1.42 fm , 
the binding energy per nucleon is 15.75 MeV and M + g s cf) = 0.56M. The last number 
corresponds to a huge attraction that is nearly canceled by the huge repulsion. Then one 
may use Eq. ( |5.6| ) to obtain the separate contributions of the vector mesons and nucleons, 
with spectacular results. The use of Eq. ( |5.6j ) leads immediately to the unusual result that 
nucleons carry only 65% of the plus-momentum. Thus is much less than the 90% needed to 



explain the EMC effect for infinite nuclear matter ||45|| . According to Eq. (|5.6|) the vector 
mesons must carry the remaining 35% of the plus-momentum, an amazingly large number. 

It is instructive to evaluate the vector Uy = gy V~ and scalar Us potentials that the 
nucleon feels with the present parameters of QHD1. These are given by 

U v = %p B = 330 MeV (5.10) 

Us = — ^tPs = —420 MeV. (5.11) 

m s 

These very large potentials, obtained with my = 783 MeV and ms = 550 MeV, are the 



distinctive features of the Walecka model The vector and scalar meson fields are given 
according to Eqs. ( |5.1| ) and fl5T2] ) as 

V° = Uy/gy = 28.3 MeV (5.12) 
= —43.9 MeV. (5.13) 

The nucleonic momentum distribution is the input to calculations of the nuclear structure 
functions. This distribution function can be computed from the integrand of Eq.(|5.6|). The 
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probability that a nucleoli has plus momentum k + is determined from the condition that 
the plus momentum carried by nucleons, P^, is given by P^/A = J dk + k + f(k + ), where 
A = Pb&- It is convenient to use the dimensionless variable y = ~ with M = M — 15.75 
MeV. Then we find the result: 



3 M 3 



M 2 V M y > 



(5.14) 



where y ± = and E F = ^k 2 F + (M + g s <f)) 2 . This function is displayed in Fig. @. The 

average value of y,(y) can be computed from this distribution: 

(y) = J™dyyf(y) = ^. (5.15) 

The relation to experiments is obtained by recalling that the nuclear structure function 
F 2 a can be obtained from the light front distribution function f(y) (which gives the proba- 
bility for a nucleon to have a plus momentum fraction y) and the nucleon structure function 
F 2 n using the relation: 



F 2A (x) 
A 



dyf(y)F 2N (x/y), (5.16) 



where y is A times the fraction of the nuclear plus-momentum carried by the nucleon, and 
x is the Bjorken variable computed using the nuclear mass minus the binding energy. This 
formula is the expression of the usual convolution model, with validity determined by a 
number of assumptions. Eq. (|5.16| ) is essentially Eq. (5.2) of Frankfurt and Strikman, |§, 
with a correspondence between our A f(y) and their / d 2 p±p^(a,p±)/a. The use of the light 
front formalism enables us to calculate the function f(y) from the integrand of Eq.( [5.6|) . 

The distribution of the vector meson plus-momentum is also an interesting quantity. The 
mean fields 0, are constants in space and time. Thus V~ has support only for k + = 0. 
The physical interpretation of this is that an infinite number of mesons carry a vanishingly 
small e of the plus-momentum, but the product is 35%. One can also show |23| that 



k + f v (k+) = 0.35M5(k + ). (5.17) 

There is an important phenomenological consequence the value k + = corresponds to 
x Bj = which can not be reached in experiments. This means one can't use the momentum 
sum rule as a phenomenological tool to analyze deep inelastic scattering data to determine 
the different contributions to the plus-momentum. 

Of course these noteworthy results are caused by solving a simple model for a simple 
system with a simple mean field approximation. It is necessary to ask if any of the qualitative 
features of the present results will persist in more detailed treatments. A simple first step is 
to include the effects of scalar meson self-coupling. We examine the parameters of Ref. 



and find one set which seems consistent with deep inelastic scattering data. This is n = 2500 
MeV, A =200, ms = 368.6 MeV. Here kp = 1.3 fm _1 and the binding energy per nucleon 
is again 15.75 MeV. The other parameters are g 2 = 27.96, g 2 = 52.44. The value of 
M + g s (p = 0.84M. This corresponds to nucleons carrying 90% of the plus momentum, and 
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FIG. 2. Light front momentum distribution as a function of y for linear and non-linear models. 



non— linear 



(k + ) = 0.10M<5(V 



(5.18) 



which is roughly consistent with the data [J5|. The nucleon momentum distribution is also 
shown in Fig. 0. The vector and scalar fields are now given by 



V 



9v 

—PB 



rn 



v 



97 MeV 
-29MeV. 



(5.19) 
(5.20) 



VI. LIGHT FRONT STATIC SOURCE MODEL FOR NUCLEAR MESONS 



The previously discussed calculations of Ref. [ 22|1 , using a Lagrangian in which Dirac 



nucleons are coupled to massive scalar and vector mesons |H2[ , treated the example of infinite 



nuclear matter within the mean field approximation. In this case, the scalar and vector 
meson fields are constants in both space and time, but only the vector mesons carry a 
non-zero plus-momentum. The constant nature means that the momentum distribution has 
support only at k + = 0. Such a distribution would not be accessible experimentally, so that 
the suppression of the plus-momentum of valence quarks would not imply the existence of 
a corresponding testable enhancement of anti-quarks. However, it is necessary to ask if the 
result is only a artifact of the infinite nuclear size and of the mean field approximation. 

In section, I review a previous investigation [23[] of the dependence on nuclear size. The 
technique was to construct the scalar and vector meson states for a model in which the 
nucleus is represented as a static source of radius R and mass Ma- In this case, Ma is 
very large for large A, and the nucleus acts as an extended static source. In the mean field 
approximation, the nucleus acts as a static source if it is infinitely heavy. Then a coherent 
state is the ground state of the Hamiltonian. 
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A. Scalar Meson Distribution 



We take the scalar mesons, as coupled to a large static nucleus represented by a scalar 
source J(r). Such a system is described by the Lagrangian density 



rn 



(6.1) 



To be specific, consider a spherical nucleus of radius R with constant density, so that 

J(r) = JoQ(R-M). (6.2) 

The quantity Jo can be thought of as arising from the product of a scalar meson coupling 
constant gs and an appropriate scalar density ps with 



Jo = 9sPs- 



(6.3) 



This is a static source in the rest frame, expressed in the usual (t, r) coordinates. 

The essential observation for the LF formulation of systems coupled to static sources is 
that static sources in a rest-frame correspond to uniformly moving sources in a LF framework 
p7| . In particular, a static charge distribution in a rest-frame Jrf(x±,x 3 ) corresponds to a 
uniformly moving charge distribution on the LF, where at x + = 



J LF (x±,x ) = J RF (x ±J z = -x /2). 
Using the usual technique to construct T Ml/ , one finds 



k+,k A _ 

where 
jRF{k±, k + ) 



k+ 



1 



d 2 x± 



k+k 



dx 



V2k+ 



jRF(k±, k + )a ] k+] *^ + h.c. 



(6.4) 



(6.5) 



- ikj_ -x± 



(2tt) 3 / 2 

1 4tt # 3 3ji{kR) 



g s ps®(R- Jx 2 ± + (x-/2 



(2tt) 3 / 2 3 



kR 



with 



k= \ k\ +k+ z 



(6.6) 
(6.7) 

(6.8) 



The interaction depends only on a combination of momenta that corresponds to the mag- 
nitude of a vector k with a positive z component, k + . The form factor 3jl fc ^ fi - > can be well 
approximated [||] by the expression e ~ k2R2 / 10 so that the typical momenta are bounded 
by ^/TO/R. The term proportional to k + arises from the plus momentum of the heavy 
source and its inclusion is necessary to maintain the rotational invariance of the light-front 
approach, see Sect. III. 
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The LF Hamiltonian P [Eq. (|6.5| )1 is quadratic in the scalar fields and its ground state 
is a coherent state 



\i>0)LF OC 



n ex p 



( 2J* RF (k ± ,k+)al +% 



k+,k i 



m 2 +k 2 . 



+ k+ 



\ 



|0>. 



(6.9) 



The action of Eq. ( |6.5|) on \iI>q)lf yields again the same state with eigenvalue 



Pn 



dk + 
o ~~k+ 



d 2 k. 



\J R F(k±,k+)\< 



m^+kj 
k^~ 



k+ 



(6.10) 



We replace the integral J °° dk + by | /f^, ciA^ so that 



Pa 



■(gsps) 



d 3 k ( 4ir R 3 3ji(kR) 



1 



(2tt) 3 V 3 
Integration of the above yields the result that 



kR 



,4tt R 3 



rrii + k 2 



(6.11) 



lim^oo P = 
The LF-momentum distribution for the scalar mesons 



(9sps) -= — 2 • 
3 m§ 



Ps(&±,fc + ) = (V>o I a k+l : ± a 
can be calculated using Eq. (|6]9|), with the result 



k+k i 



(6.12) 



(6.13) 



p s (k ± ,k^ 



2k" 



J RF (k ± ,k + ) 



[m| + P]' 



(6.14) 



Since JRF{k± ,k + ) is strongly peaked for k ~ the momentum distribution is also trivially 
peaked near small momenta (for i? 3> 1/ms). However, vanishes at fc + = 0, so that it 
necessarily is very small. Note that the light frame momentum distribution function and the 
related plus-momentum distribution can only be obtained using the light front formulation. 

The plus-momentum carried by the meson field is of great interest here. We compute 
this using Eq. ( |6.14|) , to obtain 



(k + ) = \jd 3 k 



k' 2 



jRF{k) 



771% + k 2 



(6.15) 



where we used rotational invariance of the source in the rest-frame. It is interesting to 
express the quantity (k + ) in terms of a coordinate space integral. We find: 



= d 3 rd 3 r'J(r) 



1 



e -m s \r-r'\ 

Yl-nW — f'\ 



J(f'). 



(6.16) 
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The quantity in brackets has a volume integral of (with the integration variable as \r — f'\). 
Thus the integral receives non-vanishing contributions only from regions near the nuclear 
surface, as is expected from the notion that the scalar meson field would be constant for a 
nucleus of infinite size and so receives nonzero Fourier components only from the regions near 
the surface of the nucleus. This means that (k + ) oc R 2 which gives a far smaller magnitude 
than the R 3 behavior of the binding energy. In particular, 

-rY ~ R' (6 7) 



which vanishes in the limit R — > oo, in accord with Ref. |22fl which shows that (k + ) vanishes 
for the case of infinite nuclear matter. 

A detailed investigation of the scalar meson momentum distribution,/?^), probability 
that a scalar meson carries a momentum fraction 

x = w N ( 6 - 18) 



is given in Ref. [23 . 



B. Vector Meson Distribution 



The calculation of vector meson distributions is based on the formalism in Ref. 22 



which used the light front quantization procedure of Ref. P^fTOj] . The model we consider is 
defined by taking the vector mesons to be coupled to a large nuclear source of baryon 
current J M . Thus the relevant Lagrangian density Ly is given by 



(6.19) 



where of Eq. ( ^.33| ) enters in the interaction term, with V + 
Note that 



jit 



LF 



where 



(6.20) 



J RF 



g v p B e(R-r), 



(6.21) 



because there is no special direction in space. The nucleus is in its rest frame. 

The main difference between vector mesons and scalar mesons is the appearance of the 



polarization vector e M in the coupling of vector mesons to the nucleon current (|4.34 ). The 
transverse components of the nucleon current vanish in our model, and e + = 0, so that the 
only non-zero component of e M is \i = — . 



The coherent state takes the form 



n ex p 

k+.k, 



^J^k^k+Y^aik,^ 



J 



|0>, 



(6.22) 
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where 



RF 



9vPb 




Ji(kR) 
kR 



(6.23) 



with k 



k+\ The expressions for the energy and momentum distribution are ob- 



tained by taking expectation values using this wave function. The expressions obtained are 
similar to those for the scalar mesons except that a factor of (e - ) 2 summed over all polar- 
ization states is present. The polarization sum is given by Eq. (|4.35|), with // = — , and 



k' 



k+ 



so that 



x)=l,3 



(k, u)e (fc, u) = 4- 



m 2 v + k\ 



k+< 



(6.24) 



We shall see that the denominator factor k +2 plays a large role in causing the vector meson 
momentum distribution function to be large for small values of k + . 

The light front Hamiltonian P~ consists of three terms: the kinetic energy (Eq. (2.24) of 
Ref. p2|); the linear coupling and, the effects of the instantaneous vector meson exchange. 
Thus the light-front energy of the nucleus due to coupling to the vector meson field consists 
of a contribution which arises from physical vector meson intermediate states, and an in- 
stantaneous interaction which exactly cancels the most infrared singular terms of the term 
due to dynamical mesons. This term is obtained by canonical light front quantization and 
e.g. is included in Eq. (2.48) of Ref. P2"| . One takes the matrix element of the light front 
Hamiltonian to obtain the ground state energy of the vector meson field coupled to the fixed 
source: 



Po(V) 



°°dk^ 



k + 



d 2 k, 



jRF(k_l 



k+) 



k+ 



+ k+ 



(6.25) 



which is related to the vector meson contribution to the binding energy Pq , which except 
for the sign (reflecting the fact that scalar meson give rise to attraction, while vector mesons 
give rise to repulsion between nucleons) is of the same form as the scalar result Eq. ( |6.10|) . 

The result of taking matrix elements of the plus-momentum operator in the coherent 
state leads to the result: 



p v (k ± ,k^ 



J 



RF 



2k + 



k + niy + k 2 



[m 2 v + +k 2 \ 



(6.26) 



The second term on the right hand side of Eq. ( |6.26| ) is (except for the sign and a 
differing mass) identical to the momentum distribution of scalar mesons. The contribution 
of this second term to the momentum (per nucleon) carried by the vector mesons vanishes 
in the nuclear matter limit and explicit numerical evaluation shows that for finite nuclei it 
is negligible compared with the the first term. This singular term 



J RF 



m v + k 2 



(6.27) 
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is more interesting since it diverges as k + — ► 0. By direct comparison one can verify that the 
contribution from this term to the total momentum carried by the vector mesons is identical 
to the contribution of the vector mesons to the rest-frame energy of the nucleus 

2 

roc r rco r 2 J/jf 

/ dk + d 2 k ±Pv ng k + = / dk+ d 2 k± 9 1 =P-(V). (6.28) 
Jo J Jo J mir + k l 



10 J jo J m v 

This means that in the infinite nuclear matter limit the momentum carried by the vector 
mesons is given by 

(k v ) = P (V). (6.29) 

This result that the momentum carried by vector mesons equals the P~ due to vector mesons 
holds regardless whether or not there is also a scalar interaction present. Both scalar and 
vector interaction are rather large in nuclei (of opposite sign, such that their net effect is 
small), so vector mesons may carry a substantial fraction of the nucleus' momentum. In 
the previous section, we have shown that scalar mesons carry only a small fraction so the 
net plus-momentum carried by the mesons is essentially the potentially very large plus- 
momentum carried by the vector mesons. 

The result that vector mesons carry a large part of the nuclear plus-momentum, which 
was obtained in Ref. |22[ using more general arguments, is at first surprising since the 
vector meson field (very much like the scalar meson field) is constant in space and time for 
nuclear matter in the mean field approximation. Therefore, one would expect that the vector 
meson field for nuclear matter contains only quanta with vanishing + momentum. It thus 
seems paradoxical that vector mesons nevertheless carry a finite fraction of the nucleus' 
+ momentum in this limit. To resolve this apparent paradox, we study the momentum 
distribution arising from the crucial singular piece p v ng of Eq. (|6.27|) . 

A +- momentum distribution function per nucleon (A = ^i? 3 p B ))is defined using the 
variable of Eq. (|6.18|) . Then, the use of Eq. (|6.27|) leads to the result 



Mx) = / d 2 k ±Pv ng (k ± , k + ), (6.30) 



,362 f°° dy 



or 

xf v (x) = R A -g v p B I — 2 ~ 2m j't(v)- ( 6 - 31 ) 
7r JxM n r y y z + m v R z 

A qualitative understanding may be gained by noting that the quantity j 2 (y)/y peaks at 
y « 1.6, and that my R and R are very large numbers. Thus, in the limit R — > oo the 
integral vanishes unless x = 0. In the limit of infinite R there is a sharp distinction between 
the results for x = and for non-zero values, no matter how small. If x = 0, the integral 
can be done. The net result is that 

hm xf v (x) = (^) 2 {f6(x), (6.32) 

in accord with the result expected from earlier work. The present meaning of the function 
S(x) is that integrals over x including this delta function are non- vanishing provided the 
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0.1 

X 



FIG. 3. The quantity xfy{x) vs. x for nuclei of radii 3,4,5,6, 10 fm. 



lower limit is infinitesimally close to zero. We follow Ref. \TQ and again use the parameters 



of of Chin and Walecka | 4"jfl . Then the use of Eq. (|6.32| ) leads to (xy) = 0.348, which agrees 
with the result of Ref. 



The distribution function fv(%) of Eq. (|6.31|) for finite values of R is shown in Fig 



These results show the typical behavior of distributions that approach delta functions. The 
approach to the delta function limit is slow. All of the results show a significant spread and 
there is support for values of x such that x greater than about 0.1 for < R < 6 fm. The 
key physics question to be addressed is whether or not the distribution is non-zero for values 
of x that are too small to be observable. This is discussed in Sect. 11. 

The scalar and vector meson states for a mean field model of large nuclei have been 
constructed using coherent states, and the meson distribution functions have been obtained. 
The light-front momentum distribution for scalar mesons in the mean field approximation 
is localized near k + — > 0, and the total momentum carried by the scalar mesons in the mean 
field approximation is vanishingly small for nuclear matter. This calculation thus confirms 
the results of Ref. P^fl : even though scalar mesons contribute to the nuclear binding in 
nuclear matter, they carry only a vanishing fraction of the momentum in the mean field 
calculation. 



Also in accord with Ref. [^] is the present result that vector mesons do contribute to the 
plus-momentum of the nucleus in the same limit of infinite nuclear matter. The momentum 
fraction carried by vector mesons in mean field approximation and in the nuclear matter 
limit is given by the ratio between the the vector meson contribution to the potential energy 
and the nuclear mass. This result is obtained even though the vector meson distribution 
functions are zero for non-zero x for infinitely large nuclei. This is because xfv{x) oc 5(x). 

For nuclear radii R corresponding to realistic large nuclei, the vector meson distributions, 
while strongly peaked at low values of x, are wide enough so that some fraction of the vector 
mesons could be observable. Thus the limitation of the support to k + = found in Ref. |22 



does not occur. However, the predictive power of the momentum sum rule is still vitiated 
because a significant fraction of the mesons are hidden at small values of xbj- 
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VII. MEAN FIELD THEORY FOR FINITE-SIZED NUCLEI 



It is important to make calculations for finite nuclei because all laboratory experiments 
are done for such targets or projectiles. The most remarkable feature of all of nuclear physics 
is that the shell model is able to explain the magic numbers. Rotational invariance causes 
the 2j + 1 degeneracy of the single particle orbitals, and full occupation leads to increased 
binding. But light front dynamics does not make rotational invariance manifest because 
the different components x~ , x± of the spatial variable are treated differently. There is only 
an explicit invariance for rotations of the xy plane about the z axis. However, any set of 
correct final results must respect rotational invariance. The challenge of making successful 
calculations of the properties of finite nuclei is important to us. 



In the following we review the work of Ref. |24[ Let's start with a schematic discussion 
of how it is that we will be able to find spectra which have the correct number of degenerate 
states. The concern here is with systems with a large number of nucleons, so that applying 
the usual light front Fock state expansion is too difficult. Instead, a variational calculation, 
a light front version of Hartree-Fock theory, is needed. The light front Hamiltonian, is P~, 
but minimizing the expectation value of that operator would lead to nonsensical results 
because P~ = M\j 'P+ . One can reach zero energy by letting P + be infinite. This difficulty 
is avoided by performing a constrained variation, in which the total light front momentum 
is fixed by including a Lagrange multiplier term proportional to the total momentum in the 
light front Hamiltonian. We minimize the expectation value of P + subject to the condition 
that the expectation values of P~ and P + are equal. This is the same as minimizing the 
expectation value of the average of P~ and P + . 

We use a simple example to demonstrate the need to include the plus-momentum along 
with the minus momentum. Consider a nucleus of A nucleons of momentum P^ = Ma, 
Pal = 0, which consists of a nucleon of momentum (p + ,p±), and a residual (A—l) nucleon 
system which must have momentum (Pj — p + , —p±). The kinetic energy K is given by the 
expression 

p\ + M 2 p\ + M\_ x 

K = ^^ + Pt-p+ ■ (7 - 1} 

In the second expression, one is tempted to neglect the term p + in comparison with Pj[ ph 
M4. This would be a mistake. Instead expand 

-^^+P + + M A _ 1 , (7.2) 
pi- 
using M\_ x jP\ pa 1 for large A. The second term is the plus- momentum operator mentioned 
above. For free particles, of ordinary three momentum p one has E 2 (p) = p 2 + m 2 and 
p + = E(p) + p 3 , so that 



( E 2 (p) - (p 3 ^ 
E(p)+p 3 



K « V ^\ , V V ' + E{p) + p 3 + M A -i = 2E(p) + M A -i, (7.3) 
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We see that K depends only on the magnitude of a three-momentum and rotational invari- 
ance is restored. The physical mechanism of this restoration is the inclusion of the recoil 
kinetic energy of the residual nucleus. 



A. Nucleon Mode Equation 



The key feature of Ref. P4J was the derivation of the equation that defines the nucleon 
single-particle modes. This was done by minimizing the value of P~ subject to the constraint 
that the expectation values of P + and P~ were equal. These equations are given in Eq. (4) 
of the short paper, and (3.28) of the long paper. Here we present a more intuitive derivation 



21] of the nucleon single particle wave equation. The first step is to develop the light front 



formalism for potentials that are static in the nuclear rest frame. Let us start from the Dirac 
equation for a static potential in the usual (t, r) (ET) coordinates 



i^dp -m- g s V s if) - gvl°V (r) 4> 



\ET 



0. 



(7.4) 



Since 7 = (7 + + 7 )/2 and since couplings using the "bad" component 7 are difficult to 
handle in the light front framework, we perform the transformation 



if, 



ET 



where 



<9 3 r = v , 

and T does not depend on time. Using the transformation ( |7. 5| ) , one finds 
= [z 7 ^ + *7± • 0±T) -m- g s V s {r) - g v (7 + 7 3 ) V Q [f)] tf, 
which may be expressed in terms of light front coordinates (x± = r_j_) as 
'1 



x 



x 



-V s (x ± 



x + — X 



1 + V (x ± 



X 



X 



2 
0. 



m 



(7.5) 



(7.6) 



(7.7) 



(7.8) 



Even though the potential is static in the equal time formulation, the Dirac equation for 
the same potential in light-front coordinates is light front - "time", i.e. x + , dependent. This 
is because a static source in a rest-frame corresponds to a uniformly moving source on the 
light-front. Given that the time dependence of the external fields is only due to a uniform 
translation, we transform Eq. ( [7.8|) into a form which contains only static (with respect to 
x + ) potentials. For this purpose, we consider the equation of motion satisfied by Dirac fields 
which are obtained by an x + (light front -time) dependent translation 



ip'(x±,x ,x + ) = e tx+p+/2 tp(x±,x ,x + ) 



(7.9) 
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The techniques of Sect. Ill, in particular Eq. ( |3.6|) , are used to show that the equation 
of motion for ip takes the form: 



1 1 x~ 

-ry + (<9 - d + ) + -17 d + + i7 ± ■ (d± + ig v T(x±, — — )) - m 

-9sV s {x ± , -^-) - 9vl + V~{x ± , = 0, 



(7.10) 



with V~ = Vq. The translated fields satisfy an equation of motion with potentials that do 
not depend on x + . Moreover, the static potentials evaluated at r correspond to light front 
potentials evaluated at (x±, — An even simpler derivation of this can be obtained from 
evaluating z = (x + — x~)/2 at x + = which says z = —x~/2. 

That the result (|7. 10| ) is the same as the equations for ip± in Ref. f24|. can be seen by 
making a decomposition into a dynamical and a constraint equation. Multiplication of Eq. 
( [7.101 ) by 7 + from the left yields a constraint equation 



id+^j_ = \ia ± ■ ± + ig v (d±T)) +(3m + V s ] ^ 



(7.11) 



where if)± = |7°7 ± '?/ ; - Multiplication of Eq. (|7.10|) by 7 from the left yields an equation 
for ip + : 



i{8- - d + - ig v V-)if>+ = ia ± ■ (d ± + igv(d ± V)) +(3m + V s V 



(7.12) 



One may use the constraint equation ( 7.11 ) to eliminate ?/>_ in Eq. ( 7.12| ) to obtain the equa- 
tion of motion for for the dynamical degrees of freedoms. The results ( [7.11 ) and ( 7.12 ) are 
the desired equations. The difference between these and the original fermion field equations 
( |4.32j ) is the appearance of the term — id + which appears on the left side of the equation, 



B. Dynamical Meson Fields 

The light-front Schroedinger equation for the complete nuclear ground-state wave func- 
tion I \&) is 

p- I = M A I (7.13) 

We choose to work in the nuclear rest frame so that we also need 

P+ I *) = M A I (7.14) 

As explained above, one must minimize the expectation value of P~ subject to the condition 
that the expectation value of P + is equal to the expectation value of P~ . This is the same as 
minimizing the average of P~ and P + , which is the rest-frame energy of the entire system. 
To this end we define a light-front Hamiltonian 

H LF = 1 -(P+ + P-). (7.15) 
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We stress that Hlf is not usual the Hamiltonian, because the light-front quantization is 
used to define all of the operators that enter. 

The wave function | 1 J r ) consists of a Slater determinant of nucleon fields | $) times a 
mesonic portion 

| *) =| $)<g> | mesons), (7.16) 
and the mean field approximation is characterized by the replacements 

(j)^ | <f> | 

V" -> (# | V | (7.17) 

The essential difference of our approach and the procedure to handle the original Walecka 
model is that the fields are treated as dynamical objects. The equations for the meson fields 
are derived in terms of expectation values of creation and destruction operators. This means 
that other expectation values involving meson fields, such as the plus-momentum density, 
also have non-zero values. 

We detail the treatment of the expectation value of <f>(x). Meson self-coupling terms are 
ignored here. Consider the quantity H LF a(k) \ ty), and use commutators to obtain 

H LF a(k) | *) = [H LF ,a(k)\ | *) + M A a(k) | tf). (7.18) 

The operators Pf of Eqs. (|4.52|) and ( 4.53|) and the standard commutation relations allow 
one to obtain 

[H LF , a(k)] = - fei + / V +m ' a(fc) + J{k \ , (7.19) 

L ' V n 2k+ V ; (27r) 3 /2 v /2^' V 7 



where 



Jik) -kPF,a((k)\, (7.20) 



(2^)3/2^21+ 2' 



and Pf is given in Eq. (|4.47|) . 

We use Eqs. ( f4.48| )-( ^.50| ), and take the commutator of the interactions Vi with a(k). 
Then re-express the results in terms of £ to obtain 

J(k) = ~\g s J d 2 x ± dx~ e ik - x Z(x)£(x). (7.21) 
Take the overlap of Eq. ( |7.18| ) with (ty \ to find 

V 1 1 ; 1 1 2k+ (2tt) 3 / 2 v / 2T+ 1 ^ 

Multiply Eq. ( |7.22| ) by a factor e~ tk x and add the result to its complex conjugate. 

The integral of the resulting equation over all k± and k + > leads to the result 
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| J <Pk±dk+0(k+) (j(k)e- lk - x + j\k)e + ^ x ) | ¥). (7.23) 



The evaluation of the right-hand-side of Eq. ( |7.23| ) proceeds by using Eq. ( |7.21| ) and its 



complex conjugate. The combination of those two terms allows one to remove the factor 
9{k + ) and obtain a delta function from the momentum integral. That \k + x~ appears in the 
exponential leads to the removal of the factor | of Eq. ( [7.21]) . One can also change variables 
using 

x = (z,x ± ). (7.24) 



2 



The minus sign enters to remove the minus sign between the two terms of the factor k ■ x. 
Then one may use a simple definition 



Note Ref. |23| obtained the relation ( |7.24j ) by examining the space-time diagram for a static 



(independent of x°) source. The net result is that 

(-V 2 + m 2 \ | <f>(x) | vl/) = -g a (* | £(x)£(x) | (7.26) 



which has the same form as the equation in the usual equal-time formulation. Note that the 
right hand side of Eq. (|7.26 ) should be a function of \x\ for the spherical nuclei of our present 



concern. Our formalism for the nucleon fields uses x± and x~ as independent variables, so 
that obtaining numerically scalar and vector nucleon densities that depend only x\ + (x~ /2) 2 
will provide a central, vital test of our procedures and mean field theory. This does occur 
24j, so the scalar field (\1> | <p{x) | depends only |ac| according to ( |7.26| ). 

We stress that the use of Eq. ( |7.24| ) is merely a convenient way to simplify the calculation 
- using it allows us to treat the _L and minus spatial variables on the same footing, and 
to maintain explicit rotational invariance. We will obtain the mesonic plus-momentum 
distributions from the ground state expectation value of operators expressed in terms of 
light front coordinates. 

The procedure of Eqs. ( |7.18j ) to ( |7.26j ) can also be applied to the vector fields, P^| . After 



some manipulations mandated by the appearance of the barred vector potential, one finds 
(-V 2 + m 2 v ) (* | V^x) | *> = g v (y | C{x)l^(x) \ *). (7.27) 



C. Finite Nucleus Solutions, Results and interpretation of the eigenvalues 

The main content of the nucleon mode equation is that it is the usual Dirac equation 
expressed in terms of light front variables. It should not be a surprise that the mode equation 
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of the ET theory turns out to be the same as that of the light front theory. In the former, 
one minimizes the energy of the chosen single determinant wave function. In the light front 
one minimizes P~ subject to the constraint that the expectation value of P + is the same as 
that of P~ . For the correct wave function, each of P ± is one half the mass or energy of the 
nucleus. 

The nuclear wave function | consists of a Slater determinant of nucleon fields | <£>) times 
a mesonic portion, and the mean field approximation is characterized by the replacements: 
4> — ► | 4> I ^ r } ) V^ 1 — > I V^ 1 I We quantize the nucleon fields using 



ip(x) = J2(x ,x ± I n)e iPnX+/2 b n 



[7.28) 



in which the variable x represents both the spatial variables x±,x~ and the spin, isospin 
indices, and b n obeys the usual anti-commutation relations. The meson fields are also treated 
as functions of these variables. The Slater determinant | $) is defined by allowing A nucleon 
states to be occupied. 

The use of Eq. ( ^28|) in Eqs. ( fTTTD and (|7l2|) leads to 



Pn I n)+ 
id + | n)„ 
I n )x 



(id + + 2g v V~) | n)+ + (a± ■ (pj 
(ax ■ (Px -9vV ± )+ (3{M + g s q 
A± | n), 



9vV_ 



+ (3{M + g s <j)))\n)„ 



(7.29) 

(7.30) 
(7.31) 



which together with Eqs. (|7.26|) and (|7.27|) forms the self-consistent set of equations to be 
solved. The first step of the solution procedure is to use a representation |49| in which 
| n)± are each represented as Pauli spinors. The light front quantization respects manifest 
rotational invariance for rotations about the z-axis. Thus each single-particle state has a 
good J z . We use a momentum representation for the longitudinal variable in which the 
values of p + take on those of a discrete set: p+ = (2m + l)7r/(2L) where m > and L is 
a quantization length. This means that the nucleon wave functions have support only for 
p + > 0. This spectrum condition is a requirement for exact solutions for any theory. The 
coordinate space representation is used for the _L variables. Then we have 



"x ± )e 



m 1 



-1/2)0 
-1/2)0 



(7.32) 



in which the upper (lower) entrees of the Pauli spinor correspond to m s = l/2(— 1/2), and 
x x = (x±,(j)). When Eq. (|7.32|) is used in Eq. (|7.12|) , one finds that the equations do not 
depend on the magnitude of J z \ solutions for ± a given magnitude of J z are degenerate. The 
functions U m (x±) and L m (x±) are expanded in a basis of B-splines of degree five in x± [|50 . 
The technical aspects of the solution procedure are detailed in Refs. [E3 , so I concentrate 



on summarizing the results. If our solutions are to have any relevance, they should respect 
rotational invariance. The success in achieving this is examined in Table I which gives our 
results for the spectra of 16 and 40 Ca. Scalar and vector meson parameters are taken 
from Horowitz and Serot [tjl|, and we have ignored electromagnetic effects. Note that the 
J z = ±1/2 spectrum contains the eigenvalues of all states, since all states must have a 
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ET 



LF 



State n p~/2 - M N (MeV) J z = ±1/2 J z = ±3/2 J z = ±5/2 

Os^ -55.40 -55.39 

0p 3/2 -38.90 -38.90 -38.90 

Opi/2 -33.18 -33.18 

0d 5/2 -22.75 -22.75 -22.75 -22.74 

lsi/a -14.39 -14.36 

0d 3/2 -13.87 -13.87 -13.88 



TABLE I. Comparison of the single particle spectra of Ca in the equal time (ET) formalism 
(e n - M N ) with the light front (LF) method (p~/2 - M N ). 



J z = ±1/2 component. Furthermore, the expected degeneracies among states with different 
values of J z are reproduced numerically. This is an essential success of our procedure. 

The obtained eigenvalues of the nucleon mode equation are essentially the same as the 
single particle energies of the ET formalism, to within the expected numerical accuracy of 
our program. The origin of this remarkable feature is the large value of the product of the 
nucleon mass and the nuclear radius. This is discussed next. 

We may understand the near equality of single particle eigenvalues using an analytic 
argument, which is essentially the inverse of the one used to begin this section. First we use 
a representation in which V M appears in the equation for the nucleon field. That is, define 
(x | n)' = e~ l9vA ^(x | n), with d + A = V°. Then multiply the equation Q7.29 ) for | n)' + by 



7 + and the equation (|7.30|) for | n)'_ by 7 . Adding the resulting two equations gives 

( 7 °(p; - 2-y g v V° - 7 3 (2p + - p-/2)) \ n)' = 2(<y ± ■ p± + M + g s <P) \ n)' . (7.33) 

Solving Eq. ( |7.33[ ) in coordinate space is expected to lead to solutions in which the spectrum 
condition is not respected exactly. Equation (|7.33|) may be converted into a manifestly 
rotationally invariant equation using the relation: x~ = — 2z, so that p + = id + = 2i 



—1 



d_ 

dz ■ 



_§_ 

dx~ 



The operator p + acts as a p 3 operator, and the result ( [7.33Q looks like the Dirac 
equation of the equal time (ET) formulation (of eigenvalue p~/2) except for an offending term 
—p~/2 multiplying the 7 3 . This term may be eliminated by including a phase factor: (x \ 
n)' = e tPn Z//2 (x I n) et- The result is that (x \ n) et satisfies the standard Dirac equation of the 
equal time formulation. This means that the eigenvalues e n = p~/2 must be approximately 
the same in the two formulations. The net result of these transformations is that if one 
were to include p + < in the LF calculation then the eigenvalues e n in the LF and ET 
calculations would be identical and the LF wave function would be expressed as 

(p + ,x ± I n) = -= dze- i{p+ - p " /2)z e i9vA(x ^ z) {z,x ± | n) et- (7.34) 

\/2,7T J — 00 



The relation ( [7.34 ) tells us that {p + , x± \ n) peaks at p + pa M — gyV , with a width of the 



order of the inverse of the radius of the entire nucleus. Therefore, neglecting p + < is only a 
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Nucleon p + distribution 




y 

FIG. 4. Nucleon plus-momentum distribution function 



minor approximation, which is the reason why our above LF calculation, with p + < gives 
P„ /2 that are approximately equal to the ET results. 

Eq. ( [7.34 ) shows that the influence of the vector potential is to remove plus-momentum 
from the nucleons. Furthermore, the large value of the nuclear radius causes the region of 
support to be very narrow, so that (p + ,x± \ n) is very small for negative values of p + . 



D. Nuclear momentum content and lepton- nucleus deep inelastic scattering 



Table 7.2 |24| gives the contributions to the total P + momentum from the nucleons, 
scalar mesons, and vector mesons for 16 0, 40 Ca, and 80 Zr, as well as the nuclear matter limit. 
The vector mesons carry approximately 30% of the nuclear plus-momentum. The technical 
reason for the difference with the scalar mesons (which have negligible effect) is that the 
evaluation of a*(k,u))a(k,u) counts vector mesons "in the air" and the resulting expression 
contains polarization vectors that give a factor of tr which enhances the distribution of 
vector mesons of low k + . The results for the nucleon and vector meson plus-momentum 



distributions are shown in Figs. 7.1 and 7.2 [24]. As the size of the nucleus increases the 



enhancement of the distribution at lower values of k + becomes more evident. 

It is worthwhile to see how the present results are related to lepton-nucleus deep in- 
elastic scattering experiments. We find that the nucleons carry only about 70% of the 
plus-momentum. The use of our in standard convolution formulae lead to a reduction in 
the nuclear structure function that is far too large (~95% is needed fTj) to account for the 
reduction observed J7| in the vicinity of x ~ 0.5. The reason for this is that the quantity 
M + g s (f) acts as a very small nucleon effective mass of about 670 MeV. While such a low 
value is needed to reproduce the nuclear spin orbit force, it causes difficulties, not only in 
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TABLE II. Nuclear momentum content 16 0, 40 Ca, 80 Zr, and nuclear matter (NM) in MeV. 
Contributions from nucleons, P^/A, vector mesons, P+/A, and scalar mesons Pj~ /A are displayed 
along with the sum P + /A 



Nucleus P+/A P+/A P+/A P+/A 

704.7 6A 221.8 932.9 

40 Ca 672.6 4.7 253.3 930.6 

80 Zr 655.2 3.6 270.2 929.0 

NM 569.0 0.0 354.2 923.2 
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deep inelastic scattering, but also in understanding the quasi-elastic (e, e') reaction |52|.The 
use of other Lagrangians [53-55] will lead to improved results. Including effects beyond the 
mean field lead to a significant effective tensor coupling of the iso-scalar vector meson |56f . 
and to an increased value of the effective mass. Such effects are incorporated in Brueckner 
theory, and a light-front version ]2S] could be applied to finite nuclei with better success in 
reproducing the data. 

A simple way to improve the phenomenology is to modify C, by for example including 
scalar meson self coupling terms: 3 , 4 . A wide variety of parameter sets reproduce the 
binding energy and density of nuclear nuclear matter |4"5[] . For one set, nucleons carry 90% 
of P + , so that vector mesons carry 10%, see Sect. 5. This could be acceptable. There is a 
problem with this parameter set, the related nuclear spin-orbit splitting is found to be too 
small ||57|| . This is not so bad, since there are a variety of non-mean field mechanisms which 
can supply a spin orbit force. Thus one finds a need to go beyond mean field theory. This 
involves the introduction of light front Brueckner theory. 



VIII. CHIRAL SYMMETRY AND PION-NUCLEON SCATTERING 

The essence of Brueckner theory is the use of an in-medium nucleon-nucleon scattering 
matrix to derive the self-consistent fields. This scattering matrix is derived from a realistic 
nucleon-nucleon potential. This means that the influence of pions and approximate chiral 
symmetry must be incorporated. In this section, we discuss the light -front quantization of 
a chiral Lagrangian, and review its application to low energy pion-nucleon scattering. 

The starting point is the definition of a non-linear chiral model in which the nuclear 
constituents are nucleons ip (or ?/>'), pions n scalar mesons <ft an d vector mesons V^. The 
Lagrangian without the influence of chiral symmetry is given above in Eq. ([OB). The 
additional pieces which include pions and implement chiral symmetry is C x , with 

C x = +\f 2 Tr(d^ U d» f/ f ) + -m 2 J 2 Tr(U + - 2) - $ ((M + g s cf>)(U - 1)) i//. (8.1) 
The unitary matrix U can be chosen from amongst three forms Uf 

U x =j»™lt, U 2 = \ + %l5T ' 7Z, l S U 3 = /L - TT 2 // 2 + i 75 r • TT//, (8.2) 

1 - 275T • tt/2/ v 

which correspond to different definitions of the fields. 

The pion-nucleon coupling here is chosen as that of linear representations of chiral sym- 
metry used by Gursey ||58|| , with the the Lagrangian approximately (m n 7^ 0) invariant under 
the chiral transformation 

U -> e' il5T ' a U e - i75X ' a . (8.3) 
One may transform the fermion fields, by taking U 1 ^ 2 ^' as the nucleon field. One then gets 



Lagrangians of the non-linear representation [59]. In this case the early soft pion theorems 



are manifest in the Lagrangian, and the linear pion-fermion coupling is of the pseudovector 
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type. However, the use of light front theory, requires that one find an easy way to solve 
the constraint equation that governs the fermion field. The constraint is handled in a usual 
fashion if Eq. ( |S.1| ) is used, and we'll show that this contains the early soft pion theorems. 

The constant M- plays the role of the bare pion-nucleon coupling constant. If / is chosen 
to be the pion decay constant, the Goldberger-Trieman relation yields the result that the 
axial vector coupling constant Qa = 1, which would be a problem for the Lagrangian, unless 
loop effects can make up the needed 25% effect. Corrections of that size are typical of order 
(y-) 3 effects found in the cloudy bag model |48| for many observables, including g&. 



There are no explicit A's in the above Lagrangian. We note that treating the higher 
order effects of the pion-nucleon inherent in this Lagrangian is likely to lead to a resonance 
in the (3,3) channel of pion-nucleon scattering. Such resonance effects can be included in 
the two-pion exchange contribution to nucleon-nucleon scattering. However, the A really 
is dominated by its 3-quark component [J|8| and carrying out such computations seems 
cumbersome. Thus it would be useful to incorporate the A in the Lagrangian. 

The scalar meson of Eq. (|4.26 ) is kept, even though the effects of n — n interactions, 



which could lead to similar effects, are now included in the Lagrangian. We follow many 
authors and include a scalar meson to simplify calculations and neglect the n — n interactions 
of Eq. ( J8.1| ). Note that, in the present Lagrangian, the scalar meson is not a chiral partner 
of the pion- the chiral transformation is that of Eq. ( |3.3| ). 

It is now worthwhile to mention a subtle feature regarding chiral symmetry in light front 
formalisms. Chiral invariance is defined as invariance under the transformation defined by 
Eq. ( |8.3|) if the equal time formalism is used. Now the independent fermion variable is ip+ and 
ip- is a functional of this. Thus chiral invariance is the invariance under the transformation 

V>+ - e*»™V+, (8-4) 
which is not the same as Eq.(|8.3|) pUlpl. However, the T + ~ (or light front Hamiltonian 



density) derived from £ + C x is invariant under the transformation (|8.4j ) if the pion mass is 
neglected [22]|, so the usual chiral properties are obtained in these light front dynamics. 



We begin by showing that, if one starts with a non-linear representation of chiral sym- 
metry, the requirement of solving the constraint equation for the minus component of the 
fermion field leads one to a Lagrangian of the Gursey-type linear representation. 

The focus is on chiral properties and pion-nucleon scattering, so we momentarily dispense 
with the vector and non-chiral <fi meson fields, and examine only the following fermion-pion 
term of a non-linear representation: 



Cn« = n 



jjd» -M + * — l—^^T-d^iz - (^t) 2 7 m t • TZXd^TZ 

1 + (vr/2/) 2 12/ ' ' K 2f J 1 



N. (8.5) 



The most general non-linear realization does not specify the ratio of the two pion-nucleon 
couplings. Next obtain the fermion field equation and make the usual decomposition: N± = 
A±N with 

(id- - Cr) N + = [ a± ■ (p ± - ± ) + (3M]N^ 
(id + - + ) iV_ = [oi • - 0±) + (3M]N + , (8.6) 
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in which the operator M is defined as 



0» 



-7) 



We need to remove the + term from the left hand side of the equation for N_. This can 
be done by defining, in analogy with the Soper-Yan transformation, a unitary operator F 
and fermion field \ such that 



with 



N = F X 



id + F = + F. 



The identity [|58 



(8.8) 



(8.9) 



5.10) 



where U 2 is given in Eq. (|8.2|), helps a good deal. Its use in Eq. (fg.9|), combined with the 
condition d^{U2U 2 ) =0, leads to the result 



U. 



2 ■ 



(8.11) 



so that using Eqs. (|8.11|) and ( p.8| ) in ( |8.6| ) yields 



id'x+ = [a± ■ p± + f3MU 2 ] x- 
id + X- = [a±-p± + PMU 2 } X+ . 



(8.12) 



This is of the desired form in which no interactions appear on the left-hand-side of the 
equation for x~- Thus the use of light front quantization mandates that the pion-nucleon 
interactions be of the form of Eq. ( |8.1|) , in which the possibility of using U\ or U3 arises from 
a field re-definition of that Lagrangian. 

The new chiral terms of C x lead to additions to the potentials vi,v 2 given above in 
Eqs. (|4.48| , |4.4g| ). These additions are: 



; x = J d 2 x ± dx-£(U - 1)(M + g s <f>)£, 



(8.13) 



and 



d 2 x±dx £ 
J d 2 x±dx~ £ 



{U-l)(M + g 



r 



- 7 

9vl-V 



2p^ 



2p+ 

(U-l)(M + g, 



g v7 -V + (U-l)(M + g t 



(8.14) 



The term v\ accounts for the emission or absorption of any number of pions through the 
operator U — 1. The term v 2 includes contact terms in which there is propagation of an 
instantaneous fermion. 
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The first test for any chiral formalism is to reproduce the early soft pion theorems |]6T 
Here we concentrate on low energy pion-nucleon scattering because of its relation to the 
nucleon-nucleon force. We work to second order in 1// in this first application. In this case, 
each of the [/,■ takes the same form: 



U = l + i7s- 



7T 



1 



:7T 



/ V 2 

This expression is to be used in the potentials v± and °f Eqs. 
Substitution leads to the approximations 



(8.15) 

U3D and ( gill ). 



M 
J 



J d 2 x±dx £ 



r/ 5 T • 7T - — 7T 



and 



M 2 , , , 
— / d x ± dx £ 



«7 5 T • 7T 



7" 



2p^ 



Z7 5 T • 7T 



(8.16) 



(8.17) 



The second-order scattering graphs are of three types and are shown as time x + ordered per- 
turbation theory diagrams in Fig. 8.1. The kinematics are such that 7r(q)N(k) — > 7r(q')N(k'), 
with Pi = q + k and Pf = q' + k'. The iteration of v i to second order yields the direct and 
crossed graphs of Fig. 8.1a. In this formalism v * is proportional to the matrix element of 75 
between u spinors, so it is proportional to the momentum of the absorbed or emitted pion. 
Thus the terms of Fig. 8.1a vanish near threshold. The terms of Fig. 8.1b are generated by 
the M75f terms of v±. Using the field expansions: 



7T (x) 



d 2 k ± dk + 9(k^ 



(2n)y 2 V2k^ 



a(k)e 



-ik-x 



a\k)e 



ik-x 



(8.18) 



and 



fa) 



d 2 k ± dk + 6{k^ 
(2tt) 3 / 2 v / 2^ 



2 [u(k, X)e- ik - x b{k, A) + v(k, X)e +lk - x d\k, A)] , (8.19) 



A=+, 



in the expression (|8.16|) for leads to the result that plus-momentum is conserved and the 
plus momentum of every particle is greater than zero. This means that the first of Fig. 8.1b 
vanishes identically and the second vanishes for values of the initial pion plus momentum 
that are less than twice the nucleon mass. The net result is that only the instantaneous 
term of v% and the 7r 2 term of v* (shown in Fig. 8.1c) remain to be evaluated. 
Proceeding more formally, we evaluate the S-matrix given by 



S 



rp -i P° dx+p-(x+) 

i+e 2 J -°° 1 , 



(8.20) 



where T + is the x + (light-front time) ordering operator and Pj is the interaction represen- 
tation light front Hamiltonian. Then 



-27ri5(Pr-p-)(f\T(Pr)\i), 



(8.21) 
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with 



T(P-) = Pf + Pf——-T(Pr) (8.22) 

The evaluation proceeds by using the field expansions in the expressions for v\ and v\. 
Integrating over d 2 x±dx + and evaluating the result between the relevant initial and final 
pion-nucleon states leads to the result that each contribution to the S-matrix is proportional 
to a common factor, 

2{2nf^k>+k+q>+q+' 

which combines with the result of the required integration over the light cone time to 
provide the necessary momentum conservation and flux factors. The remaining factor of 
each term is its contribution to the invariant amplitude A4. The result is 

M 2 u(k'h + u(k) M 2 u(k'h + u(k) M 

= ^'TWtM + TfTi T W^H " ^p u{k)u{k) (8 - 23) 

where the three terms here correspond to the three terms of Fig. 8.1c At threshold, k' + = 
k + = M and q' + = q + = m^, so the role of cancellations in the reduction of the term 
proportional to Sif is immediately apparent. In more detail, one finds 

, , r 2m 2 m n M 

M = 5 *f-ff + 2ie /in r n -^- (8.24) 

to leading order in m n /M. The weak nature of the 5if term and the presence of the second 
Weinberg- Tomazowa term is the hallmark of chiral symmetry 



The same results could be obtained using the linear sigma model, with a exchange playing 
the role of the 7r 2 term of Eq. (|8.15| ). 



IX. NUCLEON-NUCLEON SCATTERING ON THE LIGHT FRONT 



The correlations between nucleons are caused by the nucleon-nucleon interaction. Thus 
a necessary first step towards a light-front theory of nuclear correlations is the derivation 
of a light-front theory of the nucleon-nucleon interaction. Previous work K!3 
the light-front version of the Lippmann-Schwinger equation, the Weinberg equation, 



showed that 



can 

be transformed (except that the retardation effect is kept) into the Blankenbecler-Sugar 
equation f78fl . Kinematic invariance under boosts in the three-direction is maintained, and 
a one-boson exchange potential, is in reasonably good agreement with the NN phase shifts, 
is obtained |25 . 

how using the light-front Hamiltonian of 



It is worthwhile to begin by reviewing p2]J21 
Eqs. ( 4.4£ - 4~oH| ) and ( 8.13| ), ( |B.14| ) leads to the one-boson exchange potential. This derivation 
is useful in understanding the full nuclear wave function discussed in Sect. X. Consider the 
scattering process 1 + 2 — > 3 + 4, of Fig. |7| . 

The use of second-order perturbation theory shows that the lowest-order contribution to 
the nucleon-nucleon scattering amplitude is given by 
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FIG. 6. Low energy pion-nucleon scattering, x + -ordered graphs, (a) Second-order effects of 
the uj^u term vf. (b) Second-order effects of the u^v and vj^u terms of v*. (c) Effects of the 
instantaneous fermion propagation terms of , and of the it 2 term of The terms vf are defined 
in Eqs. QB.16US.no 



with 



(3, l\K\l, 2) = (3, 4\(vi + v*)g(Pr)( Vl + v?) +v 3 \l,2), 



<?o(i>7) = 7^ 



(9.1) 
(9.2) 



where P,~ is the negative component of the total initial or final momentum which are the 
same. In constructing the NN potential one uses conservation of four-momentum between 
the initial and final NN states in constructing the NN potential. The expression ( |9.1| ) yields 
a one-boson exchange approximation to the nucleon-nucleon potential. 

It is worthwhile to discuss the energy denominator P^ — P ~ in more detail. To be 
specific, suppose that kf > k 3 . Then the emitted meson of mass /i has momentum k with 



k + = kf — k 3 , kj_ = kix — k 3± and k 



k+ 



Pn 



P\2 Pn 



34 



. Then 



{k x k 3 ) 



k ± 2 + /x 2 

L+ _ L+ ' 



(9.3) 
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FIG. 7. x + -ordered graphs One boson exchange contributions to nucleon-nucleon 
scattering, x + -ordered graphs. The numbers 1-4 represent the momentum, spin and charge states 
of the nucleons. Here kf > k%. (a) meson propagation terms (b) instantaneous vector meson 
exchange of V3, Eq. (2.48) 



The interaction K also contains a factor of k + in the denominator, so that the relevant 
denominator is 

D = k + (P^ - P Q ~) = (kf - k+)(K ~ K) -k 2 ± -n 2 = q 2 - fi\ (9.4) 

with 

q = h- k 3 . (9.5) 

This last familiar form involves the four-momentum transfer between nucleons 1 and 3 leads 
to the Yukawa-type potentials that include the effects of retardation. It is also useful to 
explore the form of the energy denominator using light-front variables by first defining the 
plus-component, P + , of the initial and final total momentum. We may also define kf = xP + 
and k% = x'P + in which x and x' (x > x'), as ratios of plus- momenta, are invariant under 
Lorentz transformations in the three-direction. Then using ( |9.3| ), we find 
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D-^'f-^^y^-Z-S. (9.6) 

This quantity is also invariant under Lorentz transformations in the three-direction. This 
expression is to be used only if kf > k%. If k% > kf, then use a version of expression 
in which x and x' are interchanged. 

A straightforward evaluation of Eq. ( |9 . 1|) using Eqs. ( ^48 - 45^ ) leads to the result 



(3,4|#|1,2> = 2<3,4|V%2) I ' /J , (9.7) 

lb + h + h + h + 
1 2 3 4 



where 5^ 2,+ )(Pj — Pf) = 5 ( - 2 \Pi± — Pf±)6(P^~ — Pf) and V is the standard expression for 
the sum of the tt, <fi and vector meson exchange potentials: 

(3, A\V\l, 2) = (3, 4| + V(tt) + V(V)\1, 2). (9.8) 

The operator K is twice the usual two-nucleon potential times a factor which includes the 
light front phase space factor and a momentum-conserving delta function. 

For the exchange of scalar and pseudoscalar mesons, only the term (v\ +Vi)g (Pf)(vi + 
v*) enters, and one finds 

t%AWU Mio\ «(4)r«(2) u(3)Tu(l) 

(3,41^^)11,2) = AM2{2n)3 {q2 _ ^ ■ (9-9) 

The notation is that u(i) is the Dirac-spinor for a free nucleon of quantum numbers i, and 
T is either of the form g s or % g^^T. The derivation of the contribution of vector meson ex- 
change proceeds by including the meson exchange v ig (P^)vi plus the meson instantaneous 
term t> 3 , and the result takes the familiar form: 



The expressions (|9.9j) and ( |9.10| ) represent the usual |fT2]-[76]] expressions for the chosen one- 



boson exchange potentials, if no form factor effects are included. The sum of the amplitudes 
arising from each of the individual one-boson exchange terms gives the invariant amplitude 
to second order in each of the coupling constants. The factors in Eqs. (|9.9|) and (|9.10|) 

can be thought of as re- normalizing the spinors so that uu = 1, and the factors \Jjj^ of 
Eq. ( |9.7| ) serve to further change the normalization to u*u = 1 

These amplitudes are strong, so computing the nucleon-nucleon scattering amplitude 
and phase shifts requires including higher order terms. One may include a sum which gives 
unitarity by including all iterations of the two particle irreducible scattering operator K 
through intermediate two-nucleon states. One first removes kinematic factors by defining a 
T-Matrix T using 

M 2 5^(Pi- PA 
M = 2T v (9.11) 

'h + h + h + h + 
t\j<2 ' L, 3 ' L 4 
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to find that 



(3,4|T|1,2) = <3,4|V|1,2)+ £ /(3, 4jF|5, 6) fff ° i,^** . (5, 6|T|1, 2). 

(9.12) 



One realizes that Eq. ( p. 12 ) is of the form an equation, similar that of Weinberg 
derived for the scattering of nucleons by Frankfurt and Strikman p|, by expressing the 
plus-momentum variable in terms of a light-front momentum fraction a such that 

Pt = OiPt, (9-13) 
and using the relative and total momentum variables: 

k± = (1 - a)k 5± - ak 6± , 

P i± = k 5± + k 6± . (9.14) 

Then, 

(3,4|T|1,2> = (3,4^1,2)+ / £(3,4|F|5,6) ff (5,6|T|1,2), 

(9.15) 

where Pf is the square of the total initial four-momentum, otherwise known as the invariant 

energy s and ^?jz^) is the corresponding quantity for the intermediate state. Because the 
kernel V is itself invariant under Lorentz transformations in the three-direction and the 
integral involves k± and a the procedure of solving this equation gives T with the same 
invariance. Note that we use the labels ki to designate momenta in the intermediate state, 
and ki for the initial and final states. 

Equation ( 9.15 ) can in turn be re-expressed (in the center of mass frame) as the 
Blankenbecler-Sugar (BbS) equation |78] by using the variable transformation [[79]] : 

E(k) + k 3 . 

a =^4r- (9 - 16) 

with E(k) = Vk-k + M 2 , and Pf = 4(fe • k + M 2 ) The result is: 

(3,4|T|1,2) = <3,4|V|1,2> + J £ (3, 4|V|5, 6) --— _ ± ' (5, 6|T|1, 2), (9.17) 



^5,-^6 



which is the desired equation. Rotational invariance is manifestly obeyed. The three- 
dimensional propagator is exactly that of the BbS equation. There is, one difference between 
Eq. (|9.17|) and the standard BbS equation. Our one-boson exchange potentials depend on 
the square of the four momentum q 2 transferred when a meson is absorbed or emitted by a 
nucleon. Thus the energy difference between the initial and final on-shell nucleons is included 
and q° ^ 0. This non-zero value is a consequence of the invariance of D of Eq. ( |9.(j| ) under 
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Lorentz transformations in the three-direction. The usual derivation of the BbS equation 
from the Bethe-Salpeter equation specifies that q° = is used in the meson propagator. 
Including q° ^ instead of q° = increases the range of the potential relative to the usual 
treatment, and its consequences are explored below. One can convert Eq. (|9.17|) into the 
Lippmann-Schwinger equation of non-relativistic scattering theory by removing the factor 



M/E(p) with a simple transformation [30 



A. Realistic One-Boson Exchange Potential 



So far we have reviewed how the light front technique is used to derive nucleon-nucleon 
potentials in the one-boson exchange (OBE) approximation and use these in an appropriate 
wave equation ( |9.17| ). It is also true that this procedure yields potentials essentially identical 
to the Bonn OBEP potentials [[73|,74|j and these potentials lead to a good description of the 
NN data |2§. 

The Bonn one-boson exchange potentials employ six different mesons, namely, tt, rj, uj, p, a 
and the (isovector scalar) 5/a meson. The present formalism can account for the 7T, rj,u) 
and a in an approximately chiral invariant manner. We wish to add in couplings ipT ■ Sip 
and ipT ■ p^'jfj.ip in a chiral invariant manner. Simply adding such terms to the Lagrangian 
of Eq. ( |6.19| ) would lead to a violation of the approximate symmetry of Eq. ( |8.3| ). However, 
one can redefine the operator U so that the symmetry remains. We replace the operator 
ijj'Uip' in the Lagrangian ( |6.19| ) by ip'Uip' : 



U = e 2 f»' 



e 2/ <s Ue 2f s e 2f p 



(9.18) 



Then the new Lagrangian is invariant under the transformation 



3 175 T-a 



JJ — ¥ e -*75T-a jj e ~i-y 5 T-a 



(9.19) 



In the applications the exponential is expanded to first order in the meson fields. 

The final term to be included is the tensor a^ v q u part of the p-nucleon interaction. 
The presence of such a tensor interaction makes it difficult to write the equation for ip_ 
as ip- = l/p + ■ ■ This is a possible problem because the standard value of the ratio 
of the tensor to vector p-nucleon coupling f p /g p is 6.1, based upon Ref. j^TJ. Reproducing 
the observed values of E\ and P-wave wave phase shifts requires a large value f p /g p ; see 
Ref. ||82|| . However our Lagrangian generates such a term via vertex correction diagrams 
(which are the origin of the anomalous magnetic moment of the electron in QED). Thus 
the procedure of Ref. p5] was to simply add the necessary tensor terms to the one boson 
exchange potential. Recently Cooke 0] has provided the necessary light front quantization 
of a Lagrangian containing the tensor a^ v q v part of the p-nucleon interaction. 



This brings us to the treatment of divergent terms used in Ref. |25| . The definition of any 
effective Lagrangian requires the specification of such a procedure. Meson (m) nucleon form 
factors F m (q 2 ) were introduced which reduce the strength of the coupling for large values of 
— q 2 . This is also the procedure of Refs. |73|-[75|. In principle, calculating the higher order 
terms using the correct Lagrangian can lead to consistent calculations of these form factors. 
We use a more phenomenological approach here. 
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FIG. 8. 1 Sq phase shifts-solid light front, dashed Thompson equation. 

The net result is that the one-boson exchange treatment of the nucleon-nucleon poten- 
tial and the T-matrix resulting from its use in the BbS equation is essentially the same 
as the one-boson exchange procedure of Refs. [72-75]. The only difference is the keeping 



of the retardation effects — the square of the four-vector momentum transfer enters in our 
potentials. 



B. Results for the two-nucleon system 



Following established procedures ]73,[74], the coupling constants and cutoff masses of the 
six OBE amplitudes are varied within reasonable limits such as to reproduce the two-nucleon 
bound state (deuteron) and the two-nucleon scattering data below the inelastic threshold 
(about 300 MeV laboratory kinetic energy). 

In Table |TJ, we show the meson parameters for the Light-Front (LF) OBEP of Ref. p5| 
together with the predictions for the deuteron as well as low-energy neutron-proton (np) 
scattering. For comparison, we also give the parameters from an OBEP that was previously 
constructed and applied in the Dirac-Brueckner approach to nuclear matter Jn|[75|. The 
latter uses the Thompson formalism [S3] which is very similar to the BbS formalism. Note 
that the Thompson OBEP uses n a — 1 also for vector meson form factors, which explains 
the differences in the vector meson cutoff masses between the two OBEP. 

Phase shifts for np scattering are shown in Fig. 1 of the long paper of Ref. J25| for all 
partial waves with J < 2. A typical result is shown here in Fig. [8|. Over- all, the reproduction 
of the NN data 



SI 



5S[ by our LF OBEP is quite satisfactory and certainly as good as by 
OBEP constructed within alternative relativistic frameworks. Based upon these results, we 
applied this OBEP to the relativistic nuclear many-body problem. 
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TABLE III. Potential parameters and predictions for the deuteron and low-energy np scatter- 
ing. For the deuteron, the binding energy Bj, the D-state probability Pjj, the quadrupole moment 
Qd, and the asymptotic D-state over 5-state ratio D/S are given. Low-energy np scattering is 
parameterized in terms of a np and r np in 1 S'o and at and r t in 3 Si, where a denotes the scattering 
length and r the effective range. The nucleon mass is M = 938.919 MeV. 

Light-Front OBEP Thompson OBEP a Empirical 6 



Meson Parameters: 





m a (MeV) 


gl/^ [f/g] 


A Q (GeV) 


SS/47T [f/g] 


A a (GeV) 


al/^ [f/g] 




138.04 


14.0 


1.2 


14.6 


1.2 


13.5 - 14.6 


V 


547.5 


3 


1.5 


5 


1.5 


< 5 


p 


769 


0.9 [6.1] 


1.85 


0.95 [6.1] 


1.3 


0.6(1) [6.6 ± 1.0] 




782 


24.5 [0.0] 


1.85 


20.0 [0.0] 


1.5 


24 ± 5 ± 7 


do 


983 


2.0723 


2.0 


3.1155 


1.5 




a 


550 


8.9602 


2.0 


8.0769 


2.0 





Deuteron 



B d (MeV) 
Pd (%) 
Qd (fm 2 ) 
D/S 



2.2245 
4.53 
0.270 c 
0.0250 



2.2247 
5.10 
0.278 c 
0.0257 



2.224575(9) 

0.2860(15) 
0.0256(4) 



Low-Energy np Scattering 



a np (fm) 
r np (fm) 
at (fm) 
n (fm) 



-23.745 
2.671 
5.494 
1.856 



-23.747 
2.664 
5.475 
1.828 



-23.748(10) 
2.75(5) 
5.424(4) 
1.759(5) 



"Potential B of Brockmann and Machleidt [75] . 

b For more comprehensive information on the empirical data and references, see Table 4.1 



and 4.2 of Ref. 0. 

c Meson-exchange current contributions not included. 
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X. CORRELATED NUCLEONS IN INFINITE NUCLEAR MATTER 



I begin with an outline of the procedure and then review the formalism of Ref . . The 



first step is to derive a light front version of the nucleon-nucleon interaction. This is most 
easily done within the framework of the one boson exchange approximation. The resulting 
nucleon-nucleon potential V(NN) describes phase shifts reasonably well, as discussed in 
Sect. IX. The corresponding interaction density is V(iViV). The basic Lagrangian density 
contains a free nucleon term Cq(N), a free meson term £ ( mesons ) an d an interaction term 
£i( N, mesons) but does not contain V(iViV). Thus one adds this term and subtracts it: 

£ = C (N) - V(NN) + C m (10.1) 
£ m = d{N, mesons) + £ (mesons) + V(NN). (10.2) 

We use the term £q(N) — V(NN) to obtain a first solution | $) to the many-body problem. 
The term £ m accounts for mesonic content of Fock space, and we present |25|] a scheme to 
incorporate the effects of £ m and calculate the full wave function | Our procedure allows 
us to assess whether or not V(NN) has been chosen well. If it has, the effects of £ m can be 
treated perturbatively. 

Solving for | $) is no easy task -it demands a special non-perturbative treatment. One 
introduces a mean field Umf which acts on single nucleons. 

£ (N) - V(NN) = £ (N) - U MF + (Umf - V(NN)) . (10.3) 

The operator Umf is chosen to minimize the effects of {^>\Umf~ V(NN)\^f). There is a well- 
known procedure, called Brueckner theory, which is used to determine Umf- I n schematic 
terms: 

Umf ~ G x p, (10.4) 

in which G is a nucleon-nucleon scattering matrix, as modified by the Pauli principle, p is 
the nuclear density, and the x represents a convolution. 

The result is a rather complete theory in which the full wave function |^) includes 
the effects of both NN correlations and explicit mesons. This theory is reviewed next. 



A. Nucleonic Truncation For The Many-Body Problem 

One starts with a Lagrangian and can derive a light front Brueckner theory from first 



principles ||25|| . The nuclear wave function for the ground state of infinite nuclear matter at 



rest is defined as and we wish to solve the equation 

P-\V) = M A \V), (10.5) 

in which P~ is the light-front Hamiltonian discussed above. For a nuclear system at rest we 
must have also 

P + \y) = M A \V}. (10.6) 
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The operator P can be separated as follows: 



P- = P (N) + J (10.7) 

in which Pq (N) is the kinetic contribution to the P~ operator, giving — ^ — for the minus- 
momentum of free fermions. The operator J is the sum of the terms of Eqs. ( 4.48|) -( f4.50 ). 
CT ) and ( CT) : 



J = Vi + v 2 + v 3 + vf + V2 (10.8) 

An approximate solution of Eq. ( |10.5| ) is obtained in two stages. First the nucleons-only 
part of the Hilbert space is constructed. This involves the assumption that using a nucleon- 
nucleon interaction K accounts for the meson-nucleon dynamics. This assumption in then 
relaxed, and the formalism necessary to construct the best possible potential and to include 
meson degrees of freedom in the wave function is reviewed. Our Hamiltonian (P~~) contains 
no important terms in which the vacuum can spontaneously emit particles. (Some of the 
meson self interaction terms do involve non-zero creation from the vacuum ||38|| , but their 
influence on the physics would enter only through a medium modifications of the nucleon 
mass which are not part of standard relativistic Brueckner theory.) This simplifying feature 
of the dominant terms causes the derivations to look very similar to those of non-relativistic 
theory, even though the treatment is relativistic. 

All of the interactions in the Lagrangian are expressed in terms of the meson-nucleon 
vertex functions and contact terms represented by the operator J. The traditional approach 
is to introduce a two-nucleon potential and temporarily eliminate the meson degrees of 
freedom. This accomplished by subtracting and adding the two-nucleon potential to the 
Lagrangian. It is worthwhile to define the P-minus operator so obtained as Pq with: 

P =P (N) + K. (10.9) 

The operator K is given in Eq. ( |9.1|) . The complete P-minus operator is given by 

P- = P-+H U (10.10) 

with 

E x = J-K + P -(m) (10.11) 

where P ~(m) accounts for the contributions to P~ of mesons which do not involve inter- 
actions with the nucleons. The formal problem of choosing the best K by minimizing the 
effects of H\ = J — K is discussed below in sub-section X.E. 

The purely nucleonic part of the full wave function is defined as | and is the solution 
of the light front Schroedinger equation 



P | $) = [P (N) + K) | $) = M | $). (10.12) 

The eigenvalue problem stated above is considerably simpler than the initial one, but does 
contain the full complications of the nuclear many-body problem. We shall next review the 
derivation of the light front Brueckner Hartree-Fock Ell approximations. 
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B. Light Front Brueckner Hartree-Fock Approximation 



The interaction K that appears in Eq. ( |10.12[ ) is strong and the scattering amplitude 



is obtained, as discussed in Section 9, by solving the light front version of the Lippmann- 
Schwinger equation which treats the interaction between two nucleons to all orders in K. 
One needs to find an all-orders treatment for the ground state. This can be done, and 
one can also find the Slater determinant | 0), recall Eq. ( 10.13 ), that leads to the best 



approximation for the energy M of the full nucleonic wave function | $): 

Pmf I 0) = P '(N) | 0) = m | 0). (10.13) 

The operator P ~ (N) is expressed in terms of nucleon field operators defined by the single 
particle basis: 

k~\i) = k~\i,) (10.14) 
k _ =kil + { M + Us?+Uv (1015) 

This is a generalization of Eq. ( |5.3| ) in which self-consistent potentials, Us and Uy to be 
defined below, replace the meson field terms, gs<p and gyV~ in the equation defining the 
nucleon mode functions. Ths details of the related spinors are given in Ref. [[25]]. 

Both of the states | 0) and | $) are eigenstates of a P-minus operator, and both are 
eigenstates of the operator Pq~(N). We shall use standard techniques to derive a perturba- 
tion theory which employs the appropriate transition matrix generated by K to obtain an 
expression for the state | $) in terms of | 0). Thus we write 

| $} = | 0) + A | $) (10.16) 

with 

A = 1-|0)(0|. (10.17) 
The use of algebra leads to the result 

|fl,) - | ^ + M„-A(n-W + A-)A AA ^>- (1 °' 18) 



We can obtain a useful expression for M by acting with the operator (0 | [P Q (N) + Kj on 
the left of Eq. ( |10.18|) and using (0 | $) = 1,. Then we find 



M = {4 , I P a ~(N) + K \ (j>) + (<j> \ KA Ma _ m \ N)+K)A AK I 0>, (10.19) 
and with Eq. (|TaT9|) 

Mo - m„ = (HK | « + | KA M> _ ^ AK | (10.20) 
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Thus 

M -mo = {(f) | X | 0), (10.21) 

with 

X = * + *A— — /ArN - AX. (10.22) 

M -AP _ (iV)A v y 

The operator X is a many-body operator acting on all nucleons via the iterations of the 
two-nucleon interaction K. We shall make the independent pair approximation of including 
only pair-wise interactions. Thus we approximate 

(0 1 x | 0) « (0 1 E f m(^7) 1 0) = (0 1 r 1 0), (io.23) 

where r^j is a two-nucleon operator which is a solution of the integral equation 

r y W) = + ^-_ A p- (iV)A ^-(^7). (10.24) 

The notation i, j refers to a pair of particles. The relevant matrix element is expressed using 
the eigenstates of Eq. (p.3|) as 



(3,41^)11,2) = (3,41*11, 2}+ £ /(3, 4|*15,6) ^ *f (5,6^1,2), 

(10.25) 

in which we define 

M* = M + U S . (10.26) 

The operator Q is the two-body version of A and projects the momenta k$ and k 6 above 
the Fermi sea. The factor — 4= ==5^ appears in each of the terms of Eq. (|10.25 ), so it 

^1 ^2 ^3 ^4 

is worthwhile to define a Bruckner G-Matrix G using 

M* 2 <5 (2 ' +) fP- - P f ) 
2C (10.27) 

ft l ^2 ^3 ft 4 



To follow the steps of Sect. IX in converting Eq. (|10.25|) into one of a more familiar form, 
in which rotational invariance is manifest, one needs to know the values of 

P 1 - 2 = K + ^ J (10.28) 

which for the case of relevance here in computing the nuclear expectation value in the inde- 
pendent pair approximation, is the same as k$ + k~l . The single-particle minus- momentum 
eigenvalues are given according to the light front Eq. (|10.15| ). Our approximation is that 
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Uy is independent of orbital i. Thus this potential cancels out in computing the differ- 
ence Pf — (k§ + kg) and the energy denominator is as in the free space considerations of 
Sect. Ill, except that the mass of the nucleon is replaced by M + Us- Thus the previous 
derivation of an equivalent three-dimensional integral equation that is manifestly covariant 
and rotationally invariant proceeds as before. 

One expresses the plus-momentum variable in terms of a light-front momentum fraction 
a of Eqs. ( |9.16| ) and fl9.14| ) so that one obtains 



(3,4|G|1,2> = (3,4|V|1,2>+ / £ <3,4|V|5,6> — 2 (5, 6|G|1, 2), (10.29) 

J , i ml — a) p2 k ± +m 

A5,A 6 > r % a{l-a) 



where Pf is square of the total initial four-momentum, computed using M + Us for the 
nucleon mass. 

Equation ( |10.29| ) can in turn be re-expressed as a medium-modified Blankenbecler-Sugar 
(BbS) equation |78| by using the medium-modified version of the variable transformation 

2E* k ' K ' 



with E% = V k 2 + M* 2 . The result is: 

r /If* 2 H 3 k O 

(3,4|G|1,2) = (3,4|V|1,2) + / £ (3,4|y|5, 6)— ^4^(5, 6|C|1, 2), (10.31) 

which is the desired equation. 

The Brueckner light front Hartree-Fock (BHF) approximation is defined by taking the 
mean field to be calculated from the average G-matrix according to 

Ui = (l\ (U s + 7°^) BHF | = £ (Tm | G \l m) a . (10.32) 

m<F 

The sum over occupied orbitals l,m gives 

£[/z = 2(0|G|0). (10.33) 

KF 

We use this BHF mean- field to determine the value of mo via Eq. ( |10.13|) . Then the use of 
Eq. (|10.23|) in Eq. (|10.21|) determines the value of M . as the eigenvalue of P ~. But M is 
also the eigenvalue (or in this case the expectation value) of P + . The minimization of P ~~ 
subject to the constraint that the expectation value of P + is the value of P ~ leads to Q5.6|) : 

£ = As f ^ ( fci + (M + + ^ )2 - 2 l - £ ^ Us <^\ + (0 | r | $ (io.34) 

n (2vr) 3 Jf 1 \ k+ 2 Y 5 ^2k+ J ^ 1 1 ^ 1 ; 

^ ' f d 2 k ± dk + k + . (10.35) 



Q (2tt) 3 Jf 
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FIG. 9. Binding energy vs density of nuclear matter 



Note that the quantity k + is defined as 

k + = E* k + k\ 



10.36) 



Taking the average of equations (|10.34|) and (|10.35|) , and using the basis of Eq. (|10.15|) leads 
to our result for the BHF version of the nuclear mass: 



M o = Y, e i~ \ E (lfh\G\l m) a , 

KF 1 l,m<F 



10.37) 



in which ei is the eigenvalue of the equal time self-consistent Dirac equation for the nuclear 
modes. This is equivalent to the usual expression of Ref . |73|j75|| . 



C. Results for Energy versus Density 

The formalism of the previous section is used to calculate the energy per nucleon in 
nuclear matter as a function of density, Eq. ( |10.37|) . Our result is plotted in Fig. || by the 
solid line. The curve has a minimum value of £/A = —14.71 MeV at hp = 1.37 fm -1 , and 
gives an incompressibility of K = 180 MeV at the minimum. These predictions agree well 
with the empirical values £/A= -16 ± 1 MeV, k F = 1.35 ± 0.05 far 1 , and K = 210 ± 30 
MeV [|7|. 

To better understand our predictions, it is useful to compare with the results from 
alternative relativistic approaches. Brockmann and Machleidt |75] predict S/A = —13.6 



MeV, kp = 1.37 fm l , and K = 250 MeV at saturation, using the equal-time formalism and 
their 'Potential B'. The greatest difference occurs for the incompressibility which is predicted 
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smaller by the light front Brueckner theory implying a softer equation of state. This can be 
partially attributed to the medium effect that comes from the meson propagators in the light 
front approach and that is absent in the equal time (ET) approach. In the ET approach the 
meson propagators are 

(10.38) 



—q 2 — m 2 a — (k' — k) 2 — vn 2 a 
In nuclear matter, the free-space LF meson propagators are replaced by 



(E'* - E*) 2 - (k 1 -k) 2 -m 



2 



(10.39) 



while the ET propagators undergo no changes. The medium effect on the LF meson propa- 
gators enhances them off-shell which leads to more binding energy. This is demonstrated in 
Fig. £J] where the difference between the dotted and solid curve is generated by the medium 
effect on the meson propagators. 

There is another difference that arises from a technical issue in the solution of the 
transcendental equation for the G-matrix. We obtain new values of the mean fields 
M*(k F ) = M N + U S = 718 MeV and U v = 165 MeV. The mean field potentials obtained 
here from the G-matrix are considerably smaller than those of mean field theory in which 
the potential is used. The implications for nuclear deep inelastic scattering are discussed in 
Sec. 10.G. 

The most important medium effect in relativistic approaches to nuclear matter comes 
from the use of in- medium Dirac spinors representing the nucleons in nuclear matter ('Dirac 
effect'). This effect (and the medium effect on meson propagators) is absent in the conven- 
tional (non-relativistic) Brueckner calculation which yields the dashed curve in Fig. [|. A 
characteristic for all non-relativistic Brueckner theory calculations are the results that the 
saturation density is too high, or the binding energy is too low. 



D. Meson degrees of Freedom 



The nucleonic wave function | $) of Eq. (|10.18| ) gives us the purely nucleonic part of 



the Fock space, in which the effects of the mesons have been replaced by the two-nucleon 
interaction K. However, the full wave function is | of Eq. ( |10.5|) . We need to assess 
whether | <3>) is a good approximation to | and to determine the mesonic plus-momentum 
distributions. 

To obtain recall the relation between the full P~ operator and the one of Eq. (|10.12|) 
(Pq = Pq(N) + K) corresponding to the nucleonic wave function | $): 

P- = P + J - K + P (m). (10.40) 

Using this in Eqs. ( |10.5| ) and ( |10.12| ) yields 



*> =l «> + Ah -Va> (J -*"«>' (10 ' 41) 
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where A$ = 1- | $)($ |. 

An expression for the nuclear mass M A can be obtained by multiplying Eq. ( 10.41|) by 
($ | P~ using ($ | = 1 to obtain the result: 

M A = ($ | P- | $) + ($ | (J - ir)A $ __i__ I -A $ ( J - K) | $}. (10.42) 

For the purely nucleonic wave function | <3>) we have 

($ |p- | $) = M Q + ($ I J-K\ (10.43) 

so that 

M A = M + ($ | J - K | $) + ($ | (J - ^)A $ — 1_—A $ (J - K) | (10.44) 

The difference between M A and Mo is the expectation value of the operator O, which satisfies 
the integral equation 

O = J-K+ (J - K)A^G (M A )A^O, (10.45) 

with 

G (M A ) = A*— 1 A«. (10.46) 

The lowest relevant order of Eq. (|10.45|) is given by 

O^J-K + (J- K)A$G (M A )(J - K). 



The one-boson exchange interaction K is given by Eq. ( |9.1| ) (which now includes also the 
effects of Sect. 9.B) and we may determine if the expectation value ($ | O | <&} is reasonably 
small. If this is true, the quantity M would be a good approximation to the true eigenvalue 
of the P~ operator Ma- In the one-boson exchange approximation 

($ | O | $} « ($ | w 3 - fT | $) + ($ | J G (M A )J I $) (10.47) 

= <$ | v 3 - IT | $) + ($ | («i + «i)Go(M A )(^ + I $). (10.48) 
The use of Eq. ( B.l[ ) yields 

($ | O | $) « ($ | ( Vl + i£) (G (M A ) - g (Pr)) ( Vl + v?) | $>, (10.49) 

where the term P~j is specified in Eqs. ( |10.25|) -( |10.15| ). But within the independent pair 
approximation (in which one includes only the energy (minus-momentum) differences for a 
chosen pair of nucleons) 

G (M A )=g (Pr), (10.50) 

and the expectation value of O vanishes. 

Thus within our approximations, it is consistent to say that the exact nuclear mass M A 
is well approximated by M . This means that we have shown that it is acceptable to remove 
the explicit mesons for calculations of the nuclear mass. The simplicity of the derivation of 
this result is made possible by the dynamical simplicity of the vacuum, which is one of the 
defining features of light front field theory. 
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E. Momentum distributions 



We now indicate how to compute the + components of momentum. Look at T ++ as 
riven by Eq. ( f4.39|) The plus momentum carried by the scalar meson is given by 



P+(0) = / d 2 k ± dk + k + a\k)a(k), (10.51) 
while that of the pion is given by 

P+(7r) = J d 2 k L dk + k + a\k) ■ a(k), (10.52) 
and that of the vector meson is given by 

]T J d 2 k ± dk + k + a\k,oo)a{k,u). (10.53) 



\u) = 

w=l,3" 



The procedure f25| . is to evaluate the expectation value of the number operators in the 



wave function of Eq. ( [10.41] ). The meson destruction operator annihilates the first term 
of that equation but finds a non-zero result when acting on the second term. Here only the 
scalar term is considered. The evaluation of (|10.51[) using (|10.41|) leads to 



P+(0) = J d 2 k ± dk + k+(<S> | J f (fc) ^ 2 J(fc) | $) 

« / d 2 k ± dk + k + (<S> | jt(fc) ^_l_W (fc) | $), (10.54) 

in which the approximation is motivated by the near equality of Ma and Mq and the expec- 
tation that the impulse approximation evaluation of the meson exchange potential is valid. 
The term J(k) is defined in Eq. ( [7.201) . 



The operator jt ■ • ■ J to be evaluated has one and two body pieces. The one-body terms 
are related to a shift in the self energy of the nucleon caused by the medium. In infinite 
nuclear matter the ratio of pairs to single nucleons is infinite so that the number density is 
well approximated by the two nucleon terms of Eq. (|10.54 ). The evaluation is simplified by 



the use of Eq. (|10.50|) and noting that the relevant matrix element is the same as occurring 
in the one-boson exchange operator K except that the denominator is squared. Thus the 
momentum density n^(fe), defined by 

P+(0) = J d 2 k ± dk + k + n^(k), (10.55) 

is given as a derivative of the scalar-meson exchange contribution to the nucleon-nucleon 
potential: 

n,{k) « -2($ | dV i^ k) | $> (10.56) 



with 
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dV^Pr^k) 



OR 



'j 



4(k) 



(M A - AP-A) 



(10.57) 



in which the notation i,j specifies that only two-nucleon contributions are included. Note 
that | $ > is the correlated ground state. Note that the expectation value is taken using 
the single particle basis specified by Eq. (|10.15|) . Note that [^5|] 



d 



k+ 



dPij {Pij — Pq 



(q 2 -m 2 ) 2 ' 



(10.58) 



This means that evaluating the plus-momentum distribution for scalar mesons is the same as 
evaluating the expression for the scalar meson contribution to the nuclear potential energy, 
except that the potential V^(fc) (the dependence on P~j is suppressed) is multiplied by the 
factor —k + /(q 2 — m 2 ^). The net result is that 



n*(fc) = £ (im | n\ m v*(k) 



(-*+) 



l,m<F 



(q 2 



-n 



I m 



I m) a , 



(10.59) 



where Qi m is the Moeller scattering operator for the two-nucleon state I m. One finds a 
similar expression for the pionic terms. 

We can get the total number of each kind of meson (except the u) using a sum rule. 
Consider the schematic form of the equation for the G-matrix 



G(P-,)^V(Pr j ) + V{P l 



AE 



10.60) 



where is a schematic representation of the propagator of Eq. ( |10.25| ). Differentiating with 
respect to the denominator P^ appearing in a given meson m exchange potential yields the 
result: 



dG m (Pr 



(1 



Q .dV m (P~ 3 



dP^ v AE' dP^ y - ■ AE 



10.61) 



in which the label m refers to the type of meson. The potential V m appearing in Eq. ( |10.60| ) 
is simply the Fourier transform of the potential V m (q). Thus an examination of Eq. (|10.59|) 
and its pionic generalization that, considering the pion for example, 



N n = jd z qn^{q)= £ (I 



m 



9G*(Pr 



a,(3<F 



OP, 



I m) J 



10.62) 



■i j 



Numerical evaluation of Eq. (|10.62|) leads to the result that N n /A = 0.05. This is smaller 
than the 18% of Friman et al |[8) because we use scalar mesons instead of intermediate A 
states to provide the bulk of the attractive force. 
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F. Apparent Puzzle Resolved 



The values of Us, Uy and N n allow us to assess the deep inelastic scattering of leptons 
from our version of the ground state of nuclear matter. Using M*(kp) = 718 MeV and 
neglecting the influence of two-particle-two-hole states to approximate f{k + ) using Eq. ([5.15 ) 
shows that nucleons carry 83% (as opposed to the 65% of mean field theory P^| ) of the 
nuclear plus momentum. This represents a vast improvement in the description of nuclear 
deep inelastic scattering. The minimum value of the ratio F2A/F2N, obtained from the 
convolution formula (|5.16| ) is increased by a factor of twenty towards the data as extrapolated 
in Ref. [45f|. But this calculation provides only a lower limit of the nucleon contribution 



because of the neglect of effects of the two-particle-two hole states. We estimate that nucleons 
with momentum greater than hp would substantially increase the computed ratio F2A/F2N 
because F 2 n{x) decreases very rapidly with increasing values of x and because M* would 
increase at high momenta. 

Turn now to the experimental information about the nuclear pionic content. The Drell- 
Yan experiment on nuclear targets |12[ showed no enhancement of nuclear pions within an 
error of about 5%-10% for their heaviest target. No substantial pionic enhancement is found 
in (p,n) reactions. Understanding this result is an important challenge to the understanding 
of nuclear dynamics [IB"]. Here we have a good description of nuclear dynamics, and our 5% 
enhancement is roughly consistent [j55], within errors, with the Drell-Yan data. 



XI. COMING ATTRACTIONS 

The preceding sections have been concerned with reviewing published work devoted to 
applications of light front dynamics with hadronic Lagrangians to heavy nuclear systems. 
Only the first studies of mean field theory and Brueckner theory have been completed. Many 
more detailed elaborations of these studies are possible and will be needed. Much work in 
applying light front dynamics to nuclear physics remains to be done. The purpose of the 
present section is to briefly discuss a few new directions of our research. 



A. Quark-composite nucleons in the nucleus 

The first new direction is to consider the nucleons as composite systems of quarks. 
So far nuclear effects enter in deep inelastic scattering only through the influence of the 
nucleon distribution function f(y), recall Eq. ( |5.16 ) However the internal wave functions 



of the nucleons and therefore structure functions Fig could be modified by the presence of 
the nuclear medium. One model which includes such effects is the quark meson coupling 
model (QMC) initiated by Guichon and developed by Thomas, Saito and collaborators. 
Here the nucleons are treated as three quarks confined in a bag or under the influence of a 
confining potential ||55| . The nucleus is bound by the exchange of mesons between quarks in 
different nucleons. This model has the nice property that the scalar and vector potentials 
are much weaker than in the usual QHD theory. This means that nucleons carry a large 
(but not unity) fraction of the nuclear momentum, so that one can describe the nuclear deep 
inelastic scattering data. We ptl have been developing a light front version of the theory. In 
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this model, the light quarks move relativistically within nucleons which move relativistically 
within the nucleus. The second part can be handled by using the techniques discussed in the 
present review. To proceed further, one needs to obtain a light front description of the three- 
quark system. The binding of two particles is well-described within light front dynamics, see 
the reviews and our Sect. 3. Therefore, as a first step, we consider the nucleon as a bound 
state of a quark and a di-quark. In a nucleus the quark moves under the influence of the 
interactions with other nucleons as well as under the confining effects of the di-quark. The 
preliminary results are that the saturation properties of the system as about the same as in 
the usual QMC model. This means that we will be able to evaluate nuclear deep inelastic 
scattering with nuclear modified versions of the nucleon structure function. 



B. Few Body Problems 

Applications of light front dynamics to nuclear physics began with studies of the deuteron 
H J28|J30[| . The existence of many previous extensive applications, was one reason behind my 
decision to concentrate on heavy nuclei. However, there are many planned experiments on 
the deuteron. The Lagrangian approach of this review could be applied to the deuteron and 
also to bound states of three particles. 

We have begun by doing some toy model calculations using a version of the Wick- 
Cutkosky model in which two heavy scalar particles are bound by the exchange of a light 
scalar particle []9T|. This model had already been applied [92[ in light front calculations, with 



considerable success in understanding rotational invariance in scattering problems. The first 
step of Ref. ||91||was to study the breaking of manifest rotational invariance, present in light 
front dynamics. There can be many bound states in this model, so a test of rotational 
invariance is to reproduce the degeneracies of the eigenvalues. We found that the use of 
an effective Hamiltonian that takes into account two meson exchanges, leads to significant 
improvement in the degree to which the computed spectra have the correct degeneracies. 
This is much improved compared to results obtained when only one meson exchange is 
included in the effective Hamiltonian. The largest improvement occurs when the states are 
weakly bound. 

A more detailed study of the ground state using carefully computed one and two boson 
exchange potentials |)3] shows that the light front techniques reproduce the results of exact 
model calculations pl| , provided the binding energy of the two particle system is small 
compared to the mass of the heavy particles. 

The next step is to use our Lagrangian to describe deuteron form factors measured in 
electron and neutrino scattering. Much of the previous work has been devoted to using 
the impulse approximation to compute form factors measurable in electron- and neutrino- 
interactions with the deuteron. In the light front this means that one evaluates the matrix 
elements of the so-called good or J + , A + components of the electromagnetic and axial- vector 
currents. We studied this some time ago |95| and found that the computed axial form factor 



is very sensitive to the choice of matrix element of A + used to extract the form factor Fa(Q 2 )- 
This is because of an inherent ambiguity arising from the spin-one nature of the deuteron: 
there are two independent matrix elements, but only one operator A + is used in the impulse 
approximation. The D-state of the deuteron was responsible for most of the numerical effect. 
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Keister has studied relativistic effects in the deuteron axial form factor and confirms 
the dominant role of the D-state in observable effects beyond the non- relativistic limit. The 
present formalism, based on the use of a Lagrangian should allow us to compute the matrix 
elements of operators other than J + ,A + . In practice, this means that including two-body 
exchange currents effects of the \ir NN) Fock-state components will be necessary. 

A needed important effort would be to apply light front dynamics to the three-nucleon 
and three-quark systems. Studies of this problem are in their infancy p7]-pp[ . One can use 
the light front version of relative momentum variables and write done the light front version 
of the Fadeev equations. These should be solvable using modern high-speed computers. 



C. Possible experimental signature-HERMES Effect 

So far we have shown that light front dynamics can be applied to the physics of heavy 
nuclei. The standard good results for binding energy and density have been reproduced 
(with some improvement in the computed value of the compressibility of nuclear matter 
The salient new finding is that vector mesons carry a significant fraction of the 
nuclear plus-momentum, and that the meson fields can be treated in a quantum fashion 
that enables other matrix elements to be computed. The use of the simplest field theories 
leads to very strong meson fields, which account for an experimentally unacceptably large 
plus-momenta. This problem is ameliorated by the use of more complicated Lagrangians or 
by the use of Brueckner theory. 

Thus it seems that nuclei can be well-described using the light front approach. But 
it would be better to go further. In particular, it would be highly desirable to find spe- 
cific experimental evidence for a previously unappreciated feature of nuclei that is most 
conveniently computed within the framework of the light front approach. The HERMES 
collaboration discovery in high energy, low momentum transfer (0.5 GeV 2 < Q 2 < 2 GeV 2 
positron-nucleus scattering that the value of R = — is enhanced (by a factor as big as 5) 



in nuclei [26] may provide the experimental evidence we desire. 

The large nuclear value of R occurs as the result of two significant effects: is enhanced 
by a factor of two or more and a? is depleted by as much as 50%. Since a small value of 



R is the signature of spin 1/2 partons ||100|| it is natural to attempt an explanation of this 



unusual finding in terms of the mesonic field of nuclei. We found [ |101|| we find that it is 



possible to reproduce the qualitative features of the data using the nuclear lu and a fields 
which are responsible for binding nuclei. A mechanism in which a virtual photon strikes a 
nuclear uj meson turning it into a a meson gives a large enhancement to a^. To reproduce 
the small value of <jt it is necessary to find a mechanism which interferes destructively 
with the dominant process of 7* conversion to a vector meson which is scattered onto its 
mass shell by the target. The process of a virtual photon striking a nuclear a meson and 
converting it into a vector meson can provide the necessary destructive interference. Both of 
these mechanisms depend on previously unmeasured photon-meson coupling constants and 
form factors, so that reproducing the HERMES data is not a necessary consequence of our 
approach. However, we are able to argue that the data do not violate laws of physics or or 
disagree with other experiments. The model of Ref. | |101|| contains many specific predictions 
which should be easy to test. 
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XII. SUMMARY 



This paper begins with a discussion of the motivation for employing light front dynamics 
to compute the wave functions of nuclei, Sect. 1. The main reason is that high energy 
experiments are best analyzed in terms of the plus- momenta of the particles involved. Sect. 2 
contains a brief description of light front dynamics which stresses the importance of Eq. ( |2.3| ) 
and the use of a Lagrangian to define the total momenta P + and P~ which are the light 
front plus-momentum and r development operators. 

Since light front technology can appear cumbersome and recondite to a new reader, 
Sect. 3 is devoted to explaining how ordinary potentials |2T| are handled using light front 
technology. The influence of relativistic corrections to means square radii of heavy nuclei is 
examined and found to vary as A~ 5 ^ 3 , Eq. ( ^.32| ). Light front nuclear momentum distribu- 
tions are examined using a toy Hulthen model and the separate p + and p± dependences are 
interesting, Fig. 3.1. 

Light front quantization of hadronic Lagrangians (without pions and chiral symmetry) 
is presented in Sect. 4. An important element is the transformation (^4.31|) of Soper f39fl 
and Yan |RJ of the Fermion fields used to separate the dependent and independent fermion 
degrees of freedom. 

The formalism is applied to infinite nuclear matter under the mean field approximation, 
Sect. 5, [p^l - Standard results for saturation properties are reproduced. The special feature 
of the light front formalism the use of the plus momentum as one of the canonical variables. 
This enables a close contact with the experimental variables used to analyse deep inelastic 
scattering and any experiment in which there is one large momentum. This feature is 
exploited here in the derivation (within the mean field approximation) of the nucleonic and 
mesonic distribution functions for infinite nuclear matter, Eqs. ( |5.14|) -( [5.18j ). The vector 
mesons are shown to carry a significant fraction of the nuclear plus momentum, but only a 
zero plus-momentum, and therefore do not participate in nuclear deep inelastic scattering 
or Drell-Yan experiments. 

This restriction to zero plus momentum is investigated in Sect. 6, using a simple static 
model p3j. The feature that a static source in the usual coordinates corresponds to a source 



moving with a constant velocity in light front coordinates is exploited. The dependence of 
the vector meson distribution function on the nuclear radius is studied, with the key result 
reproduced in Fig. 6.1. The plus momenta are no longer restricted to zero, except in the 
limiting case of infinite nuclear matter. For finite nuclei, of radius Ra the plus-momentum 
distributions are spread, with k + ~ 1/-Ra- Thus experimental verification or disproof is 
possible. 

Sect. 7 is devoted to reviewing the application of light front dynamics to finite nuclei 



24 1 . The necessary technique is to minimize expectation value of the sum P + P + . This 
leads to a new set of coupled equations ( |7.30| ) and ( [7.291 ) for the single nucleon modes. 
These depend on the meson fields of Eqs. ( |7.26|) and (|7.27|) . The most qualitatively startling 
feature emerging from the derivation is that the meson field equations ( |7.26| ) and ( [7.27] ) are 
the same as that of the usual theory, except that z of the equal-time theory translates to 
—x~/2 of the light-front version. This can be understood in a simple manner by noting that 
light-front quantization The general argument is that this result emerges from the feature 
that a static source in the usual coordinates corresponds to a source moving with a constant 
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velocity in light front coordinates, as in Sect. 6. 

Even though the meson field equations of the light-front and equal-time theories are 
the same, there are substantial and significant differences between the two theories. In our 
treatment, the mesonic fields are treated as quantum field operators. The mean field approx- 
imation is developed by replacing these operators by their expectation values in the complete 
ground state nuclear wave function. This means that the ground state wave function con- 
tains Fock terms with mesonic degrees of freedom. We can therefore compute expectation 
values other than that of the field and obtain the mesonic momentum distributions, Table 7.2 
and Fig. 7.2. This feature has been absent in standard approaches. 

We obtain an approximate solution (|7.34 ) of our nucleon mode equation. Our nucleon 
mode functions are approximately a phase factor times the usual equal-time mode functions 
(evaluated at x~ = —2z). This shows that the energy eigenvalues of the two theories should 
have very similar values. But the wave functions are different- the presence of the phase 
factor explicitly shows that the nucleons give up substantial amounts of plus momentum to 
the vector mesons. 

A new numerical technique, is introduced in Ref. |23| to solve the coupled nucleon and 
meson field equations. Our results, Table 7.1, display the expected 2j a + 1 degeneracy of 
the single nucleons levels, and the resulting binding energies are essentially the same as for 
the usual equal-time formulation. This indicates that the approximation ( 7.34 ) is valid. 

The present results for finite nuclei, which use the original Walecka model, are not 
consistent with experiments on lepton-nucleus deep inelastic scattering and (e, e') reactions. 
This is because, in 40 Ca for example, the nucleons carry only 72% of the plus momentum. 
This is a result of the quantity M + g s (j), which acts as a nucleon effective mass, is very 
small, about 670 MeV. The use of a small effective mass and a large vector potential enables 
a simple reproduction of the nuclear spin orbit force P21BTJ. However, the use of a different 



Lagrangian, including non-linear couplings between scalar mesons should provide significant 
improvement as indicated by the light front results for infinite nuclear matter, Eq. ( p.!8| ). 
Another interesting possibility would be to obtain a light-front version of the quark-meson 



coupling model |54]J55|J90|] , in which confined quarks interact by exchanging mesons with 
quarks in other nucleons. This model, also has smaller magnitudes of the scalar and vector 
potentials. 

In any case, these kinds of nuclear physics calculations can be done in a manner in which 
modern nuclear dynamics is respected, boost invariance in the z-direction is preserved, and 
in which the rotational invariance so necessary to understanding the basic features of nuclei 
is maintained. 

Sect. 8 reviews how the light front quantization of a chiral Lagrangian can be accom- 
plished. The resulting formalism can be applied to many problems of interest to nuclear 
physics. Soft pion theorems for pion-nucleon scattering, Eq. (|8.24|) are reproduced. 

Sect. 9 reviews the our light front studies [[^,[^| of nucleon-nucleon scattering. The T- 
matrix is shown to be manifestly covariant, using the one boson exchange approximation to 
the nucleon-nucleon potential. An essential feature is the the meson propagators of Eqs. ( |9.4| ) 
and ( |9.4|) contain the influence of retardation. Light front quantization is used to obtain a 
nucleon-nucleon potential which yields phase shifts in good agreement with data. See, for 
example, Fig. 9.2. 

Sect. 10 reviews our derivation PS| of the Brueckner Hartree-Fock equations. Applying 
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our light front OBEP, the nuclear matter saturation properties are reasonably well repro- 
duced. The binding energy per nucleon is 14.7 MeV and kp = 1.37fm _1 . A value of the 
compressibility, 180 MeV, that is smaller than that of alternative relativistic approaches is 
obtained. This is largely due to the light front requirement that retardation terms be kept. 
The derivation starts with a meson-baryon Lagrangian, so we are able to show that replac- 
ing the meson degrees of freedom by a NN interaction is a consistent approximation, and 
the formalism allows one to calculate corrections to this approximation in a well-organized 
manner. The simplicity of the vacuum in our light front approach is an important feature 
in allowing the derivations to proceed. The mesonic Fock space components of the nuclear 
wave function are obtained also, and aspects of the meson and nucleon plus-momentum 
distribution functions are computed. We find that there are about 0.05 excess pions per 
nucleon. 

There are promising implications of the Brueckner theory results for studies of nuclear 
deep inelastic scattering and Drell-Yan reactions. The nucleons probably carry about 90 % 
of the nuclear plus momentum, enough for rough accord [|5| with data. Furthermore having 
0.05 excess pions is consistent fl55P , within experimental errors, with the Drell-Yan data. 

A few future directions are discussed in Sect. 11. Despite the extreme length of this 
review, much work in applying light front dynamics to nuclear physics remains to be done. 
For example, only the first studies of mean field theory and Brueckner theory have been 
completed. Many more detailed elaborations for heavy nuclei, few-body systems and com- 
posite nucleons are possible and will be needed. I hope that this review will stimulate the 
reader to future research using light front dynamics. 



ACKNOWLEDGMENTS 



This review is based on work performed in collaboration with P.G. Blunden, J.R. Cooke, 
M. Burkardt, and R. Machleidt, and D.R. Phillips. This work is partially supported by the 
USDOE. I thank S.J. Brodsky, L. Frankfurt and M. Strikman for introducing me to this 
subject. My work light front dynamics began during a sabbatical hosted by the national 
INT, the SLAC theory group and the CSSM in Adelaide. 



70 



REFERENCES 



[1] P. A. M. Dirac, Rev. Mod. Phys. 21 (1949) 392. 

[2] L.L. Frankfurt and M.I. Strikman, Phys. Rep. 76, (1981) 215. 

[3] S.J. Brodsky and G.P Lepage, Phys. Rev. D22 (1982) 2157. 

[4] S. J. Brodsky, H.C. Pauli, S.S. Pinsky, Phys. Rep. 301, 299-486 (1998); Theory of 
hadrons and light-front QCD, edited by S.D. Glazek, (World Scientific, Singapore, 
1994). 

M. Burkardt Adv. Nucl. Phys. 23, 1 (1993). J. M. Namyslowski, Prog. Part. Nucl. Phys. 
14 (1985) 49; S. J. Brodsky and G. P. Lepage in Perturbative quantum chromodynam- 
ics, A. Mueller (ed.). (World Scientific, Singapore 1989); X. D. Ji Comments Nucl. 
Part. Phys. 21, 123 (1992) W.-M. Zhang, Chinese J. Phys. 32, 717 (1994). "Theory of 
hadrons and Light-Front QCD" Edited by S.D. Glazek, World Scientific, (Singapore) 
1994; R. J. Perry, in "Hadron Physics 94: Topics on the structure and interactions 
of hadronic systems", V.E. Herscovitz et al. (eds.) World Scientific, Singapore [hep 
th/9407056| . 



[5] A. Harindranath, p.l in Light-Front Quantization and Non-Perturbative QCD", ed. 
by J. P. Vary and F. Wolz, Int. Inst, of Theoretical and Applied Physics, 1997 [hep- 
ph/96122"44] . 



[6] J. Aubert et al, Phys. Lett. 123B (1982) 275; R.G. Arnold et al, Phys. Rev. Lett. 52 
(1984) 727; A. Bodek et al, Phys. Rev. Lett. 51 (1983) 534. 

[7] R.L. Jaffe, in Relativistic Dynamics and Quark-Nuclear Physics, pp. 537 ed. by M.B. 
Johnson and A. Picklesimer (Wiley, New York, 1985); L.L. Frankfurt and M.I. Strik- 
man, Phys. Rep. 160(1988) 235; M. Arneodo, Phys. Rep. 240 (1994) 301; D.F. 
Geesaman, K. Saito, A.W. Thomas, Ann. Rev. Nucl. Part. Sci. 45 (1995) 337. 

[8] L.L. Frankfurt and M.I. Strikman, Phys. Rep. 160 (1988) 235. 

[9] C.H. Llewellyn Smith B128 (1983) 107. 
[10] M. Ericson and A.W. Thomas, Phys. Lett. B128 (1983) 112. 

[11] R.P. Bickerstaff, M.C. Birse, and G.A. Miller, Phys. Rev. Lett. 53, (1984) 2532; M. 

Ericson and A.W. Thomas, Phys. Lett. 148B (1984) 191. 
[12] D.M. Aide et al, Phys. Rev. Lett. 64 (1990) 2479. 
[13] G.F. Bertsch, L. Frankfurt, and M. Strikman, Science 259 (1993) 773. 
[14] S. J. Brodsky and H. C. Pauli, SLAC-PUB-5558, Invited lectures given at 30th Schlad- 

ming Winter School in Particle Physics: Field Theory, Schladming, Austria, Feb 27 - 

Mar 8, 1991; S. J. Brodsky, G. McCartor, H. C. Pauli and S. S. Pinsky, Part. World 3 

(1993) 109. 

[15] H. C. Pauli, Eur. Phys. J. C7 (1999) 289. 

[16] K. G. Wilson, T. S. Walhout, A. Harindranath, W. Zhang, R. J. Perry and S. D. Glazek, 
Phys. Rev. D49 (1994) 6720; M. M. Brisudova, R. J. Perry and K. G. Wilson, Phys. 
Rev. Lett. 78 (1997) 1227; 

[17] P. Grange, A. . Neveu, H. C. Pauli, S. . Pinsky and E. . Werner, "New non-perturbative 
methods and quantization on the light cone" Proceedings, School, Les Houches, France, 
February 24-March 7, 1997," Berlin, Germany: Springer (1998) 309 p. 

[18] A. Casher, Phys. Rev. D14 (1976) 452. 



71 



[19] N. C. Schoonderwoerd and B. L. Bakker, Phys. Rev. D57 (1998) 4965; 

N. C. Schoonderwoerd and B. L. Bakker, Phys. Rev. D58 (1998) 025013. 
[20] S. A. Paston, V. A. Franke and E. V. Prokhvatilov, Theor. Math. Phys. 120 (1999) 

1164. 

[21] PG. Blunden, M. Burkardt and G.A. Miller, Phys. Rev. C61 (2000) 025206 
[22] G.A. Miller, Phys. Rev. C 56 (1997) R8; 56 (1997) 2789. 
[23] M. Burkardt and G.A. Miller, Phys. Rev. C 58 (1998) 2450. 

[24] PG. Blunden, M. Burkardt, and G.A. Miller, Phys. Rev. C59 (1999) R2998; Phys. 

Rev. C 60 (1999) 55211. 
[25] G.A. Miller and R. Machleidt, Phys. Lett. B455 (1999) 19; Phys. Rev. C60 (1999) 

035202-1. 

[26] "Nuclear Effects on R=ctl/o"t in Deep Inelastic Scattering", K. Ackerstaff et al, Phys. 
Lett. B (in press), |iep-ex/991007l| and DESY-99-150. 

[27] "Relativistic Nuclear Physics in the Light Front Formalism," V.R. Garsevanishvili and 
Z.R. Menteshashvili, Nova Science Publishers, Commack (1993). 

[28] F. Coester, Helv. Phys. Acta. 38, 7 (1965); and in "Constraint's Theory and Relativis- 
tic Dynamics, edited by G. Longhi and L.Lusanna (World Scientific Singapore, 1987), 
p. 159. 

[29] B. D. Keister and W. N. Polyzou Relativistic Hamiltonian dynamics in nuclear and 
particle physics in Advances in Nuclear Physics. J.W. Negele and E. Vogt (ed.). 20 
(1991) 1 (Plenum Press, New York). 

[30] V.A. Karmanov and A.V. Smirnov, Nucl. Phys. A575 (1994) 691; J. Carbonell, V.A. 
Karmanov, B. Desplanques, V.A. Karmanov J.F. Mathiot Nucl. Phys. A589 (1995) 
697; J. Carbonell, B. Desplanques, V. A. Karmanov and J. F. Mathiot, Phys. Rept. 
300 (1998) 215. 

[31] M. G. Fuda and Y. Zhang, Phys.Rev. C51 (1995) 23; Phys. Rev. C54 (1996) 495. 
[32] TJNAF HAPPEX collaboration, R. Michaels and P.A. Souder spokesmen; R. Michaels, 

talk at HAPPEX collaboration meeting, TJNAF, Dec. 16, 1998; C.J. Horowitz, 

S. J. Pollock, P. A. Souder and R. Michaels, |nucl-th/991203l . 
[33] S. D. Drell and T. Yan, Annals Phys. 66 (1971) 578. 
[34] Y. Mardor et al, Phys. Lett. B437 (1998) 257. frucl-ex/971000% 
[35] S.-J. Chang, R.G. Root, and T.-M. Yan, Phys. Rev. D 7 (1973) 1133; 7 (1973) 1147; 
[36] H.Leutwyler, J.R. Klauder, and L. Streit, Nuo. Cim. 66a (1970) 536. 
[37] F. Rohrlich, Acta. Phys. Austr. Supp. VIII (1971) 277. 
[38] M. Burkardt, Advances Nucl. Phys. 23, 1 (1996). 

[39] D.E. Soper, "Field Theories in the Infinite Momentum Frame", SLAC pub-137 (1971); 
D.E. Soper, Phys. Rev D4 (1971) 1620; J.B. Kogut and D.E. Soper, Phys. Rev. Dl 
(1971) 2901. 

[40] T-M Yan Phys. Rev. D7 (1974) 1760, 1780. 

[41] D. Mustaki, S. Pinsky, J. Shigemitsu and K. Wilson, Phys. Rev. D43 (1991) 3411. 
[42] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16 (1986) 1; Int. J. Mod. Phys. E6 
(1997) 515. 

[43] St. Glazek, C. M. Shakin, Phys. Rev. C44 (1991) 1012. 
[44] S.A. Chin and J.D. Walecka, Phys. Lett. 52B (1974) 24. 
[45] I. Sick, and D. Day, Phys. Lett. B274 (1992) 16. 



72 



[46] R.J. Furnstahl, C.E. Price, and G.E. Walker Phys Rev C36 (1987) 2590. 

[47] M. Burkardt, in "Theory of Hadrons and Light-Front QCD", Ed. S. Glazek, (World 

Scientific, Singapore, 1994). 
[48] S. Theberge, A. W. Thomas and G. A. Miller, Phys. Rev. D22 (1980) 2838; D23 (1981) 

2106(E). A. W. Thomas, S. Theberge, and G. A. Miller, Phys. Rev. D24 (1981) 216; 

A. W. Thomas, Adv. Nucl. Phys. 13 (1984) 1; G. A. Miller, p 190 in Int. Rev. Nucl. 
Phys. 2 (1984) 190. 

[49] A. Harindranath and W-M Zhang, Phys. Rev. D48 (1993) 4868, 4881, 4903. 
[50] J.A. McNeil et al, Phys. Rev. C40 (1989) 399; W.R. Johnson et al, Phys. Rev. A 37 
(1988) 307. 

[51] C.J. Horowitz and B.D. Serot, Nucl. Phys. A368 (1981) 503. 

[52] H. Kim, C.J. Horowitz, and M.R. Frank, Phys.Rev. C 51, (1995) 792. 

[53] J. Zimanyi, S.A. Moszkowski, Phys. Rev. C 42 (1990) 1416; N.K. Glendenning, F. 

Weber, and S.A. Moszkowski, Phys. Rev. C 45 (1992) 844. 
[54] P. A. Guichon, Phys. Lett. B200 (1988) 235; K. Saito and A.W. Thomas, Phys. Lett. 

327B(1994) 9; K. Saito and A. W. Thomas, Nucl. Phys. A574 (1994) 659. X. Jin and 

B. K. Jennings, Phys. Rev. C 54 (1996) 1427. K. Saito, K. Tsushima and A. W. Thomas, 
^iucl-th/99TJTU8l ; 

[55] PG. Blunden and G.A. Miller, Phys. Rev. C 54 (1996) 359. 

[56] R.J. Furnstahl, J.J. Rusnak, and B.D. Serot, Nucl. Phys. A632 (1998) 607. 

[57] R.J. Furnstahl, J.J. Rusnak and B.D. Serot, Nucl. Phys. A632 (1998) 607. 

[58] F. Gursey, Nuovo Cimento 16 (1960) 230; Ann. Phys. (N.Y.)12(1961) 91; P. Chang 

and F. Gursey, Phys. Rev. 164 (1967) 1752. 
[59] S. Weinberg, Phys. Rev. 166, 1568 (1968). 
[60] D. Mustaki, [hep-ph/ 94042UTj 

[61] B.H. Brandsden and R.G. Moorhouse, "The Pion-Nucleon System" (Princeton Univer- 
sity Press, Princeton, 1973); S.L. Adler and R.F. Dashen, "Current Algebras", (W.A. 
Benjamin, New York, 1968). 

[62] R. Blankenbecler and R. Sugar, Phys.Rev. 142 (1966) 1051. 

[63] R. Machleidt, K. Holinde, Ch. Elster, Phys. Rep. 149 (1987) 1 ; R. Machleidt, Adv. 
Nucl. Phys. 19 (1989) 189. 

[64] R. Brockmann and R. Machleidt, "The Dirac-Brueckner Approach", prepared for 
"Open Problems in Nuclear Matter", M. Baldo, ed (World Scientific, Singapore), 
Chapt. 3. |nucl-th/ 96 12001 . 

[65] Nucleon-Nucleon Interaction", G.E. Brown and A.D. Jackson North- Holland (Amster- 
dam) 1976; "Pions and Nuclei", T. Ericson and W. Weise Oxford Science Publications, 
(New York) 1987. 

[66] S. Weinberg Phys. Rev. 150 (1966) 1313. 

[67] M.V. Terent'ev, Sov. J. Nucl. Phys. 24 (1976) 106. 

[68] M-H Partovi, and E.L. Lomon, Phys. Rev. D2 (1970) 1999. 

[69] G. Hoehler and H. Pietarinen, Nucl. Phys. B95 (1975) 210. 

[70] G.E. Brown and R. Machleidt, Phys. Rev. C50 (1994) 1731. 

[71] J.R. Cooke, private communication, and J. R. Cooke and G. A. Miller, "Ground 



states of the Wick-Cutkosky model using light-front dynamics", NT@UW-99-52 [nucl- 
th/00020l6| . prepared for Phys. Rev. C. 



73 



[72] R. Machleidt, K. Holinde, Ch. Elster, Phys. Rep. 149 (1987), 1. 
[73] R. Machleidt, Adv. Nucl. Phys. 19 (1989) 189. 

[74] R. Machleidt, One-Boson-Exchange Potentials and Nucleon-Nucleon Scattering, in: 
Computational Nuclear Physics 2 — Nuclear Reactions, K. Langanke, J. A. Maruhn, 
and S. E. Koonin, eds. (Springer, New York, 1993), Chapter 1, p. 1. 
[75] R. Brockmann and R. Machleidt, Phys. Lett. 149B(1984) 283; Phys. Rev. C 42 (1990) 
1965; "The Dirac-Brueckner Approach", in "Open Problems in Nuclear Matter", M. 
Baldo, ed (World Scientific, Singapore), Chapter 3; |nucl-th/9612004 . 



[76] "Nucleon-Nucleon Interaction", G.E. Brown and A.D. Jackson, North- Holland (Ams- 
terdam) 1976; "Pions and Nuclei", T. Ericson and W. Weise, Oxford Science Publica- 
tions (New York) 1987. 

[77] S. Weinberg Phys. Rev. 150 (1966) 1313. 

[78] R. Blankenbecler and R. Sugar, Phys.Rev. 142 (1966) 1051. 

[79] M.V. Terent'ev, Sov. J. Nucl. Phys. 24 (1976) 106. 

[80] M-H Partovi, and E.L. Lomon, Phys. Rev. D2 (1970) 1999. 

[81] G. Hoehler and H. Pietarinen, Nucl. Phys. B95 (1975) 210 (1975). 

[82] G.E. Brown and R. Machleidt, Phys. Rev. C50 (1994) 1731. 

[83] R.H. Thompson, Phys. Rev. Dl (1970) 110. 

[84] V.G.J. Stoks, R.A.M. Klomp, M.C.M. Rentmeester, and J.J. de Swart, Phys. Rev. C 
48 (1993) 792. 

[85] R.A. Arndt, Virginia Polytechnic Institute (VPI) and State University, Interactive 

Dial-In Program SAID, Solution SM97. 
[86] "Structure of the Nucleus" , M.A. Preston and R.K. Bhaduri Addison- Wesley (Reading 

1975); P. Grange and M.A. Preston, Nucl. Phys.A204(1973) 1. Nucl. Phys.A219 

(1974) 266. 

[87] J.P Blaizot, Phys. Rep. 65 (1980) 171. 

[88] B. Friman, V.R. Pandharipande and R.B. Wiringa, Phys. Rev. Lett 51 (1983) 763. 
[89] H. Jung and G.A. Miller, Phys. Rev. C41(1990) 659. 
[90] G. Krein and G.A. Miller to be published 

[91] J.R. Cooke, G.A. Miller and D.R. Phillips, NTOUW-99-52; |nucl-th/99100T3| , submit- 
ted to Phys. Rev. C 

[92] N. C. Schoonderwoerd, B. L. Bakker and V. A. Karmanov, Phys. Rev. C58 (1998) 
3093. 

[93] J.R. Cooke, G.A. Miller, |riucl-th/ 0002016] to be submitted to Phys. Rev. C 
[94] T. Nieuwenhuis and J.A. Tjon, Phys. Rev. Lett. 77, 814 (1996) 
[95] T. Frederico, E.M. Henley and G.A. Miller, Nucl. Phys. A533 (1991) 617 (1991). 
[96] B.D. Keister, Phys. Rev. C55 (1997) 2171. 

[97] "The Three Body Force in the Three-Nucleon System" edited by B.L. Berman and 

B.F. Gibson (Springer New York, 1986) 
[98] L. Frankfurt and M. Strikman, in "Electromagnetic Interactions With Nuclei", B. Frois 

and I. Sick, eds. World Scientific (1991). 
[99] B.L. Bakker, L.A.Kondratyuk, M.V. Terent'ev, Nucl.Phys. B158 (1979) 497. 
[100] C.G. Callan and D.J. Gross, Phys. Rev. Lett. 22 (1969) 22. 

[101] G.A. Miller, S.J. Brodsky and M. Karliner, |hep-ph/0002T56i submitted to Phys. Lett. 
B 



74 



